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ABSTRACT 


An  adaptive  receiving  array  for  rejecting  interference  in  a  coded 
communication  system  is  described  and  computer  simulations  are  presented 
which  illustrate  basic  properties  of  the  array  processor.  The  array 
automatically  forms  a  beam  in  the  direction  of  a  signal  bearing  the 
proper  code  modulation  and  forms  nulls  in  the  directions  of  signals 
bearing  improper  modulation.  The  communication  signal  present  at 
the  array  input  is  assumed  to  be  biphase  (0°-130°)  modulated  by  the 
modulo-2  sum  of  a  periodically-repeating  code  and  a  data  sequence. 

The  code  is  assumed  known  at  the  receiving  site  and  is  the  means  by 
which  the  array  distinguishes  between  "desired"  and  "undesired"  signals. 

The  adaptive  processor  minimizes  the  mean-square  difference  (error) 
between  the  array  output  signal  and  a  reference  signal.  The  reference 
signal  contains  the  dcta  and  code  modulations  of  the  desired  input  signal 
and  is  generated  from  the  array  output  signal  by  appropriate  waveform 
processing.  The  array  weighting  coefficients  (and  hence  the  pattern)  are 
adjusted  by  a  feedback  control  system  designed  to  provide  a  steepest- 
descent  minimization  of  the  error. 

Digital  computer  simulation  results  illustrating  the  performance 
of  a  two-element  array  for  the  case  of  cw  interference  are  presented. 

The  results  show  that  signal  carrier  frequency  and  code  bit  timing 
estimates  are  required  at  the  array  and  that  the  bandwidths  of  the 
adaptive  feedback  loops  within  the  processor  must  be  restricted  for 
proper  operation.  The  simulation  results  also  show  that  the  per¬ 
formance  of  the  array  processor  degrades  gracefully  as  the  errors  in 
the  estimates  increase. 


TABLE  OF  CONTENTS 


Page 


I.  INTRODUCTION  1 

II.  THE  ADAPTIVE  ARRAY  SIGNAL  PROCESSING  TECHNIQUE  3 

A.  Basic  Adaptive  Array  Configuration  3 

B.  Reference  Signal  Requirements  5 

C.  Description  of  the  Adaptive  Array/Waveform  Processor  8 

III.  COMPUTER  SIMULATION  OF  THE  ADAPTIVE  ARRAY/WAVEFORM 

PROCESSOR  11 

A.  Simulation  of  the  Weighting  Coefficient  Loop  Equations  11 

B.  Input  Signal  Simulation  12 

C.  Simulation  of  Reference  Network  Signal  Processing  17 

D.  Sampling  Rate  Considerations  25 

E.  Reference  Signal  Phase  Quantization  26 

F.  Performance  Measure:  CW  Pattern  28 

IV.  COMPUTER  SIMULATION  RESULTS  29 

A.  Introduction  29 

B.  Effect  of  Reference  Network  Time  Delay  29 

C.  Effect  of  Code  Synchronization  Error  35 

D.  Effect  of  Desired  Signal  Frequency  Offset  39 

E.  Effect  of  Reference  Network  Processing  Gain  46 

F.  Performance  for  Large  Input  Interference  49 

G.  Results  for  Modified  Simulation  Models  58 

1.  Loop  Gain  Control  (AGC)  C8 

2.  Error  Multiplier  Output  Non-linearity  6) 

V.  SUMMARY  63 

APPENDICES  65 


REFERENCES 


81 


SECTION  I 
INTRODUCTION 


Self-phased  antenna  arrays  have  been  considered  for  use  in  a 
variety  of  applications  in  recent  years.  A  self-phased  array  coherently 
adds  the  signals  received  by  the  array  elements  regardless  of  the  angular 
position  of  the  source.  An  arbitrary  location  of  the  elements  on  any 
surface  is  acceptable  provided  shadowing  or  strong  mutual -coupling  effects 
are  avoided.  The  ability  of  the  self-phased  array  to  focus  in  the  Fresnel 
region  and  to  compensate  for  phase  disturbances  in  the  transmission 
medium  (wavefront  distortion)  are  also  significant  features. 

There  are  some  performance  limitations,  however.  For  unequal 
signal  amplitudes  at  the  element  inputs  (resulting  from  aperture  blockage 
or  mutual  coupling)  or  unequal  element  noise  powers  the  SNR  enhancement 
at  the  array  output  may  be  less  than  that  realizable  by  other  combining 
techniques. ^  The  self-phased  array  is  also  more  susceptible  to  external 
noise  or  interfering  signal  sources  than  ordinary  directive  antenna 
systems.  Under  strong  interference  conditions,  the  self-phased  array 
has  the  tendency  to  "lock'-on"  and  track  the  interference  rather  than 
the  desired  signal. 3 

In  this  report,  a  more  general  type  of  adaptive  array  processing 
is  investigated  in  which  both  the  amplitude  and  phase  of  the  element 
signals  are  adjusted  prior  to  combining.  Advantages  of  this  more 
general  approach  include  the  potential  for  diversity  combining  and  for 
interference  rejection.  Applebaum4  and  Shor5  were  among  the  first  to 
consider  adaptive  adjustment  of  element  gain  and  phase  as  an  optimal 
control  problem.  They  chose  for  an  optimization  criterion  the  maxi¬ 
mization  of  the  ratio  of  output  signal  power  to  total  noise  power. 

Widrow,  et  al6  later  analyzed  adaptive  array  systems  which  minimize 
the  mean-square  difference  (error)  between  the  array  output  signal  and 
a  reference  signal.  The  angular  position  and  modulation  components  of 
the  desired  signal  source  were  assumed  known  a  priori  at  the  receiving 
array. 

The  processing  technique  discussed  in  this  report  follows  that  of 
Widrow,  et  al  with  some  important  differences.  Knowledge  of  the  angular 
position  of  the  desired  signal  source  is  not  assumed.  In  addition,  the 
desired  signal  received  at  the  array  is  assumed  to  have  an  unknown  data 
modulation  impressed.  Information  regarding  the  desired  signal’s  carrier 
frequency,  code  modulation,  and  code  modulation  time  base  is  required  at 
the  array,  however.  The  most  significant  difference  is  the  addition  of 
a  waveform  processor  at  the  array  output.  The  purpose  of  the  waveform 
processor  is  two-fold:  1)  to  discriminate  between  a  properly  coded 
desired  signal  and  undesired  signals  in  the  array  output,  arid  2)  to 
provide  an  estimate  of  the  unknown  (data)  modulation  components.  The 
waveform  processor  generates  the  required  reference  signal  for  the 
array. 
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This  report  is  organized  as  folios:  The  adaptive  array  signal 
processing  technique  is  described  in  Sec.  II.  This  section  also  contains 
a  discussion  of  the  characteristics  required  of  the  reference  signal  and 
a  description  of  how  the  reference  signal  can  be  derived  from  the  array 
output.  In  Sec.  Ill,  the  computer  program  used  to  evaluate  the  tran¬ 
sient  performance  of  the  array  processor  is  described.  Computer 
simulation  results  are  given  in  Section  IV.  Analytical  results  which 
augment  the  simulation  results  are  included  where  they  are  available. 

The  results  are  summarized  and  conclusions  drawn  in  Section  V. 
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SECTION  II 

THE  ADAPTIVE  ARRAY  SIGNAL  PROCESSING  TECHNIQUE 


A.  Basic  Adaptive  Array  Configuration 

The  basic  configuration  of  an  N-element  adaptive  array  is  shown 
in  Fig.  1.  Here,  x-|(t),  X2\t),  xw(t)  represent  signals  at  t!e 
outputs  of  the  array  elements  and  y-j (t ) ,  y3(t),  ♦•*,  y2N-l(t)  represent 
the  processor  input  signals.  The  signal  from  each  element  is  adjusted 
in  magnitude  and  phase  prior  to  signal  combining  (summing).  The  output 
signal  from  the  kth  element  is  applied  to  two  channels:  an  "in-phase" 
channel  and  a  'quadrature"  channel.  The  quadrature  channel  signal,  y2k(0> 
is  delayed  in  phase  by  90°  with  respect  to  the  signal,  y2k-l(t)>  in  the 
in-phase  channel.  The  in-phase  and  quadrature  signals  associated  with 
the  k^  element  are  multiplied  by  weighting  coefficients  W2i<-l(t)  and 
W2k(t)i  respectively.  The  weighting  coefficients  are  real  valued  and 
can  be  positive,  negative,  or  zero.  Consequently,  the  contribution  of 
the  kth  element  to  the  array  output. 


(i)  sk(c)  =  w2k_](t)  y2k_-|(t)  +  w2k(t)y2k(t), 


can  have  an  arbitrary  amplitude  and  phase  relationship  with  respect  to 
input  signal  xj<(t).  The  array  output  signal,  S(t),  is  subtracted  from 
a  reference  signal,  R(t),  to  produce  error  signal  E(t).  The  error 
signal  and  signals  y-j(t),  •••,  y2N (t)  are  applied  as  inputs  to  the 
feedback  circuitry.  Control  voltages  generated  by  the  feedback  circuits 
determine  the  values  of  weights  w^(t). 

The  objective  of  the  processing  is  to  force  the  array  output  signal 
S(t)  to  equal  the  reference  signal  R(t).  One  treasure  of  performance  is 
the  squared  error  given  by 

?  2N  ? 

(2)  E^t)  =  [R{t)  -  l  w.(t)  y. (t)r  >  0. 

i=l  1  1 

2 

At  any  given  time,  E  is  a  quadratic  function  of  the  weights  wj ;  there¬ 
fore  E2(w] ,  •••,  W2n)  defines  a  "bowl-shaped"  surface  with  a  well- 
defined  minimum.  The  error  can  be  forced  to  assume  this  minimum  by 
appropriately  designing  the  feedback  circuitry.  One  approach  is  to  adjust 
the  weights  according  to  a  steepest-descent  minimization  procedure, 6 
i.e,,  the  wj  are  moved  in  the  maximum  downhill  direction  on  the  squared- 
error  surface.  Since  this  direction  is  given  by  the  negative  of  the 
gradient  of  E^  with  respect  to  the  w^ ,  steepest-descent  minimization 
requires  the  to  be  changed  as  follows: 


3 


(3) 


dw.  ? 

ar =  ks  VE }  ;  ks 

o 

The  larger  v7*(E  ),  the  faster  w. 
gradient  may  be  evaluated  as 


<  0. 

changes. 


The  i 


th 


component  of  the 


(4)  V.{E2)  =  2E(v.E)  =  -  2E  y., 
and  the  feedback  rule  (3)  reduces  to 

dw. (t) 

(5)  =  -2ks  E(t)  y. (t)  ; 

w.(t)  =  K  f  E(t')  y.(t')  dt '  +  w.(0) 

*  o 

where  K  =  -2k?  is  the  (positive)  integrator  gain  constant.  A  method 
for  implementing  the  feedback  circuitry  so  that  (5)  is  satisfied  is 
illustrated  in  Fig.  2  * 

B.  Reference  Signal  Requirements 

As  discussed  in  the  previous  section,  one  approach  to  optimizing 
array  performance  is  to  minimize  the  mean-square  error.  This  is  equiva¬ 
lent  to  forcing  the  array  output  signal  S(t)  to  approximate  the  reference 
signal  R(t)  as  closely  as  possible  in  a  mean-square  sense.  A  component 
of  S(t)  which  is  not  contained  in  R(t)  appears  as  an  error  signal,  and 
the  feedback  adjusts  the  weights  to  remove  it  from  the  array  output. 

The  result,  in  pattern  terminology,  is  the  formation  of  a  pattern  null 
in  the  direction  from  which  this  signal  arrives.  If  a  component  of  S(t) 
is  contained  in  R(t),  the  feedback  retains  this  signal  in  the  output  and 
adjusts  its  average  amplitude  and  phase  as  closely  as  possible  to  that 
of  R(t).  Equivalently,  a  pattern  lobe  is  formed  in  the  direction  of 
arrival  of  this  signal.  Output  errors  are  minimized  by  "accepting" 
signals  with  proper  amplitude  and  phase  modulation  and  by  rejecting 
signals  with  improper  modulation.  The  characteristics  of  the  reference 
signal,  therefore,  determine  the  ability  of  the  adaptive  array  to 
discriminate  between  "desired"  and  "undesired"  signals  at  the  array 
input.  For  proper  operation,  the-  reference  signal  waveform  must  be 
highly  correlated  with  the  desired  input  signal  waveform  and  have  a  low 
correlation  with  the  input  interference  waveform.  These  conditions  must 


*The  signal  processing  shown  in  Fig.  2  may  be  performed  at  IF  frequencies 
by  down-converting  the  received  signals  in  separate  mixers  using  a 
common  local  oscillator  signal. 
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be  satisfied  when  the  time  interval  over  which  correlation  is  measured 
is  equal  in  duration  to  the  response  time  of  the  feedback  control  loops 
and  is  chosen  arbitrarily  in  time.  These  latter  restrictions  are 
necessary  to  ensure  that  the  correlation  properties  of  the  input  signals 
are  preserved  at  the  array  output.  The  weighting  coefficients  are  time- 
varying,  in  general,  and  are  therefore  capable  of  introducing  amplitude 
and  phase  modulations  having  rates  determined  by  the  bandwidth  of  the 
feedback  loops.  The  possibility  exists  then  that  the  processor  may 
alter  the  modulation  characteristics  of  the  input  waveforms  thereby 
affecting  their  correlation  with  the  reference  signal.  The  important 
point  here  is  that  although  an  undesired  signal  may  be  completely  un¬ 
correlated  with  the  reference  signal  when  the  cross-correlation  function 
is  evaluated  over  an  infinite  time  interval,  the  response  of  the  array 
processor  is  determined  by  the  correlation  over  a  much  smaller  observation 
interval  —  in  effect,  a  moving  "window"  of  width  equal  to  the  response 
time  of  the  feedback  loops. 
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One  approach  to  satisfying  the  reference  signal  requirements  pre¬ 
sented  above  is  to  employ  wideband  signalling  techniques.  If  the  band¬ 
width  occupied  by  the  desired  input  and  reference  signals  greatly  exceeds 
■•■he  loop  processing  bandwidth,  then  these  signals  must  possess  significant 
amplitude  and/or  phase  fluctuations  over  any  time  interval  comparable  with 
the  response  time  of  the  loops.  Undesired  signals  not  containing  the 
reference  modulation  will  not  correlate  with  the  reference  signal.  The 
correlation  cannot  be  significantly  increased  by  weighting  coefficient 
modulation  in  this  case;  as  a  result,  the  feedback  processing  removes 
these  signals  from  the  array  output.  The  selection  of  a  wideband  modu¬ 
lation  technique  is  influenced  by  the  signal  environment  in  which  the 
array  operates  and  by  the  synchronization  requirement.  The  time  base  of 
the  reference  signal  modulation  must  be  aligned  with  the  time  base  of  the 
desired  input  signal  (or  vice-versa)  for  proper  correlation.  Timing 
error  results  in  a  correlation  loss  and  subsequent  partial  suppression  of 
the  desired  signal  at  the  array  output.  The  accuracy  to  which  timing  can 
be  established  and  maintained  dictates  an  upper  bound  on  the  signal 
bandwidth  since  the  cross  correlation  function  associated  with  the 
reference  and  desired  input  signals  becomes  narrower  (more  sharply  peaked) 
as  the  bandwidth  is  increased.  Other  practical  considerations  influence 
the  choice  of  a  wideband  modulation  waveform,  e.g.,  the  relative  ease  of 
waveform  generation  and  time-base  adjustment,  and  the  capability  for 
additional  processing  of  the  received  signal  to  obtain  timing  error 
information. 

When  random  (data)  modulation  is  impressed  on  the  transmitted 
signal  in  addition  to  the  deterministic  (bandspreading)  modulation, 
an  additional  problem  is  encountered  in  generating  the  reference  signal 
for  the  array.  If  the  reference  signal  does  not  contain  the  data  modu¬ 
lation,  a  loss  in  correlation  and  partial  rejection  of  the  received 
data-carrying  signal  will  necessarily  occur.  There  are  at  least  two 
approaches  to  avoiding  this  problem: 

1)  A  data-carrying  signal  and  a  deterministic  wideband  beam¬ 
steering  signal  might  be  transmitted  sequentially;  the 
operation  of  the  array  processor  would  be  correspondingly 
time-multiplexed, 

or 

2)  Techniques  for  estimating  the  unknown  (datal  modulation 
components  from  the  received  signal  itself  might  be  in¬ 
corporated  into  the  processor.  The  estimates  would  then 
be  used  in  generating  the  reference  signal. 

The  latter  approach  is  pursued  here.  An  adaptive  receiving  array  which 
employs  a  coded  reference  signal  derived  by  waveform  processing  the  array 
output  signal  is  described  in  the  next  section. 


C.  Description  of  the  Adaptive  Array/Waveform  Processor 


A  block  diagram  of  the  system  under  study  is  shown  in  Fig.  3.  The 
"desired"  signal  at  the  array  input  is  assumed  to  be  a  bi -phase  modulated, 
digital  communication  signal.  The  binary  phase  modulation  on  this  signal 
is  the  modulo-two  sum  of  a  binary  code  sequence  and  a  random  binary  data 
sequence.  The  code  is  assumed  known  at  the  array  and  is  the  means  by 
which  the  processor  distinguishes  the  "desired"  signal  from  "undesired" 
signals. 


Fig.  3.  Adaptive  array  processor  with  a 

waveform-processed  reference  signal. 


The  reference  signal  is  derived  by  waveform-processing  array  output 
S(t).  First,  a  local  oscillator  signal  is  bi-phase  .nodulated  with  the 
same  code  as  the  desired  input  signal  code  and  then  mixed  with  S(t). 

The  desired  signal  component  of  S(t)  produces  a  signal  at  the  mixer  out¬ 
put  containing  only  data  modulation  (assuming  perfect  synchronization  of 
input  and  local  oscillator  codes).  Since  the  data  rate  is  assumed  to  be 
much  less  than  the  code  rate,  this  signal  occupies  a  smaller  bandwidth 
than  S(t).  The  undesired  signals  in  S(t),  however,  do  not  undergo  this 
reduction  in  bandwidth  when  multiplied  by  the  reference  code.  Wideband 
undesired  signals  remain  wideband  at  the  mixer  output  and  narrowband  un¬ 
desired  signals  are  spread  in  bandwidth  by  the  local  oscillator  phase 
switching.  The  data  filter  passes  the  data  modulation  and  reduces  the 
power  in  the  undesired  portions  of  S(t)  relative  to  the  desired  part. 
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To  a  rough  approximation,  the  reduction  in  power  is  equal  to  the  ratio 
of  the  bandw’dths.  A  bandpass  limiter  (i.e.,  an  ideal  limiter  and 
wideband  zonal  filter)  fixes  the  reference  signal  amplitude.  Note  that 
without  limiting  the  weighting  coefficients  would  be  driven  to  zero  if 
the  reference  network  gain  were  less  than  unity  since,  in  this  case, 
the  reference  signal  magnitude  would  always  be  smaller  than  the  array 
output  signal  magnitude.  If  the  reference  network  gain  were  greater 
than  unity,  the  weighting  coefficients  would  eventually  be  forced  to 
their  saturation  limits.  The  use  of  a  limiter  circumvents  the  need 
to  maintain  unity  gain.  The  processed  signal  is  mixed  with  the  coded 
local  oscillator  signal  to  reinsert  the  code  modulation;  only  the 
sideband  centered  on  the  operating  frequency  of  the  array  appears  at 
the  output  of  the  second  mixer  (post-mixer  filtering  is  not  shown  in 
Fig.  3).  This  output,  i.e.,  the  reference  signal  R(t)  for  the  array, 
approximates  the  desired  component  in  S(t).  Slow  phase  variations  of 
the  local  oscillator  carrier  which  are  introduced  at  the  first  mixer 
are  removed  in  the  second  mixer. 

In  practice,  the  desired  compoi.ents  in  the  reference  signal  and 
the  array  output  signal  will  differ  due  to  the  delay  introduced  by  the 
bandpass  filter  and  code  synchronization  error.  Due  to  the  delay,  the 
data  transitions  are  not  synchronized.  Moreover,  the  phase  shift 
through  the  waveform  processor  is  frequency  dependent;  thus,  the  desired 
in-phase  (mod.  2tt)  condition  will  not  exist  if  the  signal  frequency  is 
offset  or  if  the  filter  parameters  change  with  time  due  to  aging  or 
temperature  variation.  Imperfect  code  timing,  i.e.,  synchronization 
error,  results  in  a  reduction  of  the  signal  to  noise  ratio  in  the  wave¬ 
form  processor,  and  the  generation  of  a  noise-like  component  in  the 
reference  signal  --  even  when  only  desired  signal  is  present  at  the  array 
output.  When  interfering  signals  are  present  in  the  array  output,  they 
are  reduced  in  strength  and  altered  considerably  in  waveform  and  fre¬ 
quency  structure  by  the  processing  operations.  However,  they  are  not 
eliminated  completely  from  the  reference  signal.  Fortunately,  complete 
suppression  of  interference  in  the  reference  signal  is  not  a  prerequisite 
for  satisfactory  interference  rejection  by  the  array.  The  essential 
requirement  is  that  the  interference  components  in  the  reference  signal 
must  be  smaller  than  the  interference  components  in  the  array  output 
signal,  i.e.,  the  effective  gain  of  the  reference  waveform  processor 
with  respect  to  interference  must  be  less  than  unity.  When  this  condition 
is  satisfied,  the  error  signal  contains  an  interference  component  equal 
to  the  difference  in  interference  levels,  and  the  feedback  processing 
responds  to  this  error  to  reduce  the  output  interference.  Reduction 
of  the  output  interference  continues  as  long  as  the  interference  component 
is  present  in  the  error  signal,  i.e.,  as  long  as  interference  is  present 
in  the  array  output  and  the  gain  of  the  waveform  processor  to  interference 
is  less  than  unity.  The  latter  condition  is  expected  to  prevail  even  for 
modest  values  of  waveform  processing  gain.  Thus,  the  array  processor  is 
potentially  capable  of  eliminating  large  interfering  signals  even  when  the 
waveform  processing  gain  is  relatively  small. 
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The  effects  of  waveform  processing  delay,  synchronization  error, 
cw  interference,  and  system  parameters  on  the  processor's  transient 
response  were  investigated  by  digital  computer  simulation.  The  simulation 
is  described  in  the  following  section. 
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SECTION  III 

COMPUTER  SIMULATION  OF  THE  ADAPTIVE  ARRAY/WAVEFORM  PROCESSOR 


A.  Simulation  of  the  Weighting  Coefficient  Loop  Equations 

A  computer  program  was  written  to  solve  the  set  of  2N  differential 
equations  (4)  for  the  weighting  coefficients.  In  sampled  form  these 
equations  are  approximated  by 


(6) 


dw.(t.)  ^  w.(t.+T)  -  w^t.) 
'  ~’dt_Ss~  =  T 


KE(V  W 


or 


(7)  w.(tj+1)  =  w.^tj)  +  KT  E(tj)  y-Ctj)  ; 

i  =  1,  2,  •••,  2N, 

where 

2N 

(8)  E(t  )  =  R(t.)  -  S(t.)  =  R(t .)  -  l  w.(t.)y.(t.) 

j  J  J  J  .j=1  I  J  *  J 

is  the  sampled  error  signal  and 

(9)  T  =  i=Vl-tj 

is  the  sampling  period  corresponding  to  a  sampling  frequency  F.  The 
time  reference  t  =  0  was  assumed  so  that 


(10)  t.  =  jT  j  =  0,1,2,.-.. 

The  yi(tj)  are  the  sampled  signals  (including  interference)  at  the  inputs 
to  the  weighting  coefficients.  The  compiler  program  performed  the 
following  sequence  of  computations  given  the  values  of  the  weighting 
coefficients  w. (t.)  at  time  t-: 

I  J  J 

a)  Computation  of  the  signal  samples  y-(t  -),  i =1 ,  *  *  *  ,2N  from 
appropriate  analytical  expressions,  J 

b)  Computation  of  the  array  output  signal  sample  S(t.)  from  the 

w^tj)  and  y.ttj),  J 
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c)  Computation  of  the  reference  signal  sample  R(tj)  from  S(tj) 
by  means  of  appropriate  digital  processing  used  to  simulate 
the  reference  generation  network. 

d)  Computation  of  the  error  signal  sample  E(t.)  from  (8),  and 

e)  Computation  of  subsequent  values  for  the  weighting  coeffi¬ 
cients  w.  (tj+-j)  from  the  recursion  equations  (7). 

The  signal  equations  y-j(t)  used  in  step  (a)  and  the  digital  processing 
used  in  step  (c)  are  discussed  in  the  following  two  sections. 

E.  Input  Signal  Simulation 

A  linear  (one-dimensional)  array  composed  of  an  even  number  N  of 
equi spaced  elements  was  assumed  in  formulating  the  equations  for  signals 
y^(t)  at  the  weighting  coefficient  inputs.  The  elements  and  weighting  coef¬ 
ficients  were  numbered  as  shown  in  Fig.  4  with  odd-numbered  elements  to 
the  right  of  the  array  phase  center  and  even-numbered  elements  to  the  left. 


X4(t)  X2(t)  Xo(t)  x,(t)  X3(t) 


WEIGHTING  COEFFICIENT  INPUTS 


Fig.  4.  Linear  array  geometry. 
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Signal  x0(t)  which  would  be  received  by  an  element  placed  at  the  array 
phase  center  was  assumed  to  be  given  by 

(11)  xQ(t)  =  sQ(t)  +  nQ(th 
with 

(12)  s 0 ( t )  =  As  C(t)  D(t)  sin  {2n(fc  +  Afg)t  +  a$0) 

(13)  nQ(t)  =  An  sin{2Tr(fc  +  Afn)t  +  an()} 
where 

As  =  desired  signal  amplitude 

An  =  interference  signal  amplitude 

f  =  array  center  frequency 

Af$  =  desired-signal  frequency  offset 

Afn  =  interference-signal  frequency  offset 

“so  =  desi red-signal  initial  (t=0)  phase 

“no  =  lnterference"s19nal  initial  phase 

C(t)  =  desired-signal  bi -phase  code  (+1,  -1  amplitude) 
modulation 

D(t)  =  desired-signal  bi-phase  data  (+1,  -1  amplitude) 
modulation . 

The  signal  received  by  the  kth  element,  x^(t),  was  expressed  as 

(14)  x^t^/  sk(t)  +  nR(t) 

s0(t  +  I  f  sin  es)  +  n0(t  +  ff  sin  6n, 

i  k  =  1,  3,  •••,  N-l 

„  so(t  -  T1  F  sin  6s>  +  "o  (t  '  IT  Z  s1n  en> 
k  =  2,  4,  •••,  N. 
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where 


-=  propagation  time  delay  between  elements 
es  =  desired-signal  arrival  angle 
=  interfering -signal  arrival  angle. 

The  propagation  time  delay  across  the  array  was  assumed  negligibly  small 
compared  to  the  code  and  data  bit  periods  of  the  desired  signal  modulation, 
i  .e . , 


(15) 

where 


N  f  «  fcd  < 


fcd  =  c*es'i red-signal  code  bit  rate 
fj  =  desired-signal  data  bit  rate. 

This  assumption  of  an  "array  bandwidth"  much  larger  than  the  input  signal 
modulation  bandwidth  justified  the  approximations 


(16) 

(17) 


C(t±ti  sin  6S)  =  C(t), 
D(t  i  jr  j  sin  es)  =  D(t). 


The  array  elements  were  assumed  to  be  ideal  (noise-free)  point  receptors.  Ac¬ 
cordingly,  the  output  signal  from  the  kth  element,  y2f<_i(t),  was  expressed  as 

(18) 


y2k-l(t)  ’ 


f 


Ac  C(t)D(t)  sin  {2ir(l  +  -r  |  f„t  +  k*  d,  (  1  +  )  sin  e  +  a  J 


(  Af 

+  An  sin  {2tiM  +  ~ 


"c 

/  Af 
frt  ♦  k,dXc  (i  ♦  ?c 


f 


so 


n 


si"  6n  +  “no1 


when 


k  =  1 ,3,5, • • • ,  N-l 
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and 


(19)  y2k_,(t)  = 


A  C(t)D(t)  sin  {2n  1  + 


Af. 


Af_  Af 

-  ft-(k-l)ird  1+^ 
Ac  Tc 


Af„ 


sin  0  +a  } 
s  so 


*•  An  sin  (2 it  1  +  ^ —  fct  -  (k-i)tf  d^  1  +  ^ —  sin  en  +  anQ} 


when  k  =  2,4,-  •  •  ,1* 


where  \ =  km  © 


is  the  element  spacing  in  wavelengths  at  the  array  center  frequency.  The 
quadrature  channel  inputs,  y2k(t),  were  assumed  to  be  given  by  expressions 
identical  to  (18)  and  (19)  except  for  the  addition  of  quadrature  phase  shifts 
to  the  sine  arguments  of  the  desired  and  interfering  signal  terms.* 

The  desired  signal  code  modulation  C(t)  consisted  of  a  length-127 
pseudo  random  sequence  generated  by  a  seven-stage  feedback  shift  register 
having  the  configuration  shown  in  Fig.  5.  The  initial  contents  of  the 
register  were  chosen  to  provide  a  representative  switching  pattern**  over 
the  first  twenty  bits  of  the  sequence. 

The  data  modulation,  D(t),  was  assumed  to  consist  of  an  alternating 
(+1 ,  -1)  sequence.  This  "square-wave"  pattern  maximizes  the  number  of  bi¬ 
phase  data  transitions  which  occur  during  the  initial  transient  period  thus 
allowing  the  effects  of  data  modulation  delay  in  the  reference  processing 
network  to  be  readily  observed. 


*The  simulation  program  actually  introduced  time  delays  of  l/4fc  seconds 
which  gave  rise  to  phase  shifts  of  90°(1  +  Af/fc)  where  Af  is  the  input 
signal  frequency  offset  from  center  frequency.  However,  the  frequency 
offsets  employed  in  the  simulation  were  small  in  comparison  with  f  so 
that  the  phase  shifts  were  approximately  equal  to  90°.  ' 

**Such  a  pattern  is  considered  to  have  approximately  equal  occurrence  of 
"l"s  and  "0"s  and  no  long  consecutive  runs  of  either  symbol. 
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CODE:  OOIO 
C(t) 


C.  Simulation  of  the  Reference  Network  Signal  Processing 


The  digital  processing  used  to  simulate  the  reference  network  is 
illustrated  in  Fig.  6.  Mixer  simulation  consisted  of  multiplication  of  the 
input  samples  by  +1  or  -1  depending  on  the  value  of  the  reference  code 
Cr(tj).  The  up/down  frequency  conversions  performed  by  the  mixers  were 
omitted  for  simplicity.  The  reference  code  was  a  delayed  version  of  the 
input  code 

(20)  Cr(tj)  =  C(t.  -  e) 

where  delay  parameter  e  could  be  varied  to  simulate  synchronization  offsets 
c  having  values  from  zero  to  one  code  bit  period  (emax’^cd  r  !)•  lerner, 
digital  bandpass  filters  were  chosen  to  simulate  the  analog  bandpass  filters 
in  the  reference  network  because  of  their  excellent  phase  .inearity  and 
reasonably  selective  passbands.  The  fourth-order  Lerner  filters  used  in 
the  simulation  are  characterized  by  the  z-transform 


(21) 


4  a.  [1  -  (e'"/2  BT  cos  b.T)  Z"1] 
HU)  =  l  1  ' 


1-1  1  -  (2e'"/2  BT  cos  b.T)  Z"1  +  e""BT  Z"2 


where  a-j  =  -  a^  =  0.5 
-a2  ■  a3  ■  1 

bl  ,4  =  2”  (fc  f  ?  1 

b2,3  *  2"<fc  »  !  > 

B  =  Bandwidth,  Hz. 

The  above  transfer  function  was  obtained  by  using  the  criteria  of  impulse 
invariance:  the  impulse  response  of  the  digital  filter  equals  the 
sampled  impulse  response  of  the  Lerner  analog  filter.  Details  of  the 
design  procedure  are  available  in  Reference  7.  Figure  7  shows  the  pole- 
zero  configurations  of  the  Lerner  filter  in  both  the  analog  (s-plane) 
and  digital  (z-plane)  domains.  The  sampled  output  of  the  filter  at  time 
t  =  nT  is  'Jven  by  a  sum  of  four  sub-outputs: 


(22)  V(nT)  =  l  a.  V. (nT), 

i=l 
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Fig.  6.  Reference  network  digital  processing. 
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Pole-zero  configuration  of  4th  order  Lerner  niter. 
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where  the  sub-outputs  are  expressed  by 


(23)  V.(nT)  =  e"ir/2  Br  cos  bj  [2V.(nT-T)  -  U(nT-T)] 

-e-irBT  V.(nT-2T)  +  U(nT). 


Function  U(nT)  represents  the  sampled  input  at  time  t  =  nT.  The  frequency 
response  of  the  Lerner  filter  was  calculated  from  (21)  with  z  =  exp(jwT) 
and  is  shown  in  Fig.  8.  Important  characteristics  to  be  noted  are  the 
ninety  degree  phase  shift  at  center  frequency  and  the  change  in  phase 
shift  over  the  bandwidth  B  of  nominally  360  degrees.  At  the  bandedge, 
i.e.,  at  f  =  fc  +  B/2,  the  amplitude  response  is  5.32  dB  smaller  than  the 
mid-band  response  and  the  phase  shift  differs  12.5  degrees  from  linear. 
Envelopes  of  the  impulse  response,  the  step  amplitude  response,  and  the 
response  to  a  180°  step  in  the  carrier  phase  are  illustrated  in  Figs.  9, 
10,  and  11,  respectively.  The  latter  figure  shows  that  the  transition 
in  the  output  signal's  phase  occurs  approximately  1/B  seconds  after  the 
phase  step  is  applied  to  the  input  signal. 


The  bandpass  limiter  in  the  reference  network  was  simulated  by  an 
ideal  limiter  characteristic  followed  by  a  Lerner  digital  filter  having  a 
passband  centered  on  the  carrier  frequency  and  a  bandwidth  B^  much  larger 
than  the  data  filter  bandwidth  Bd.  The  purpose  of  this  second  filter  was 
to  reject  frequency  components  outside  the  first  zone  from  the  output  wave¬ 
form.  The  changes  in  amplitude  experienced  by  the  signal  at  the  data  filter's 
output  following  a  180°  phase  step  due  to  the  data  modulation  (see  Fig.  11) 
are  "squared-up"  by  the  limiter  and  the  zonal  filter  introduces  a  relatively 
small  additional  delay  of  1/B^  seconds.  The  zonal  filter's  output  is  multi¬ 
plied  by  a  constant  to  result  in  a  limited  reference  signal  having  the 
desired  value:  Ar..  A  negative  constant  is  employed  for  the  scaling  factor 
so  that  the  two  ninety  degree  phase  advances  introduced  by  the  digital  filters 
(at  the  center  frequency)  are  effectively  cancelled.  This  forces  the 


desired  signal  components  at  the  waveform  processor  and  array  outputs  to 
be  in  phase  alignment  when  the  frequency  offset  equals  zero.  When  the 
frequency  offset,  f-fc,  is  sufficiently  small  so  that  the  filter  phase 
characteristics  are  approximately  linear,  the  total  phase  shift  through 
the  reference  network  is  given  approximately  by 


(24) 


Y(f):  _27r  i  +  i  {f-fc) 


-2ir(f-fc) 


B 


composite 


The  scaled  output  of  the  zonal  filter  is  multiplied  by  the  reference  code 
to  simulate  the  second  mixer.  This  completes  the  digital  simulation  of 
the  waveform  processing  network. 
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Fig.  8.  Lerner  digital  filter  frequency  response. 


D.  Sampling  Rate  Considerations 


The  selection  of  a  sampling  rate,  F,  was  influenced  by  several 
factors.  Adequate  representation  of  the  input  (bandpass)  signals  in 
sampled  form  required  that  the  sampling  rate  exceed  2fc:  the  Nyquist 
sampling  criterion.  Another  factor  to  consider  was  the  cumulative 
truncation  errors  resulting  from  the  approximate  integration  in  (7)  of 
the  weighting  coefficient  differential  equations.  These  errors  cause 
the  computer  solutions  to  differ  from  the  correct  analog  solutions.  They 
are  dependent  upon  the  rate-of-change  of  the  weighting  coefficients  and 
therefore,  from  (6),  are  related  to  input  signal  magnitude.  Approximation 
errors  are  reduced  by  increasing  the  sampling  rate.  Another  factor  arose 
from  the  limiter  nonlinearity  in  the  reference  network  which  produced  fre¬ 
quency  components  in  zones  centered  about  odd  harmonics  of  fc.  The  fre¬ 
quency  response  of  the  second  digital  filter  is  periodic,  as  illustrated 
in  Fig.  12,  with  bandpass  repetitions  centered  on  the  frequencies 

(25)  f  =  nF  ±  f 

V 

where  n  is  any  integer.  To  ensure  rejection  of  frequency  components  outside 
the  fundamental  zone  (fc),  i.e.,  to  minimize  frequency  aliasing  effects  in 
the  second  digital  filter,  the  sampling  rate  was  chosen  as  an  odd  multiple 
of  fc  as  in  Fig.  12.  The  responses  to  components  having  frequencies  around 
twice  the  sampling  frequency  (2F)  are  expected  tr  be  negligibly  small  when 
the  sampling  rate  is  sufficiently  large.  The  cot ibi nation  of  the  limiter 
non-linearity  and  a  finite  sampling  rate  also  resulted  in  reference  signal 
phase  quantization  as  will  be  discussed  in  the  next  section.  All  of  these 
factors  suggest  that  the  sampling  rate  should  be  made  much  larger  than  f  . 
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Fig.  12.  Second  digital  filter  input  spectrum  and  periodic 
frequency  response  characteristics. 


25 


In  practice,  the  sampling  rate  which  can  be  utilized  is  restricted 
by  the  length  of  time  required  to  execute  the  computer  program.  In  the 
study,  the  transient  decay  times  encountered  were  generally  on  the  order  of 
a  few  thousand  carrier  cycles  for  realistic  values  of  the  parameters. 

Program  run  time  considerations  limited  the  total  number  of  samples  processed 
to  a  few  ten  thousand.  As  a  result,  it  was  necessary  to  restrict  the  sampling 
rate  to  nine  times  the  carrier  frequency:  F  =  9f  . 

V 

E.  Reference-Signal  Phase  Quantization 

The  phase  of  the  sampled  reference  signal  was  restricted  to  a 
discrete  set  of  values  as  a  result  of  using  a  finite  sampling  rate  and 
an  ideal  limiter  characteristic.  This  phase  quantization  is  perhaps 
best  illustrated  by  an  example.  Figure  13  illustrates  the  sampled  wave¬ 
forms  at  various  points  in  the  reference  processing  network  of  Fig.  6 
when  the  array  output  signal  is  given  by 

(26)  S(t.)  =  sin(a)  tj  +  0°). 

J  c  J 

The  ninety  degree  phase  advance  in  the  data  filter  output  waveform  (b) 
results  in  a  limiter  output  waveform  (c)  composed  of  five  positive  and 
four  negative  samples  per  cycle.  The  fundamental  component  (d)  of  this 
sampled  square  wave  at  the  second  filter  output  is  changed  in  sign  and 
amplitude-scaled  in  (e)  to  form  the  sampled  reference  waveform.  Shifting 
the  phase  of  the  input  sinusoid  in  (a)  by  any  amount  up  to  plus  or  minus 
ten  degrees  changes  the  values,  but  not  the  signs,  of  the  samples  in  (b). 
Therefore,  the  sampled  waveforms  in  (c),  (d),  and  (e)  are  unchanged.  The 
reference  signal  phase  is  constant  at  zero  degrees  for  this  case.  For 
input  phase  angles  between  ten  and  thirty  degrees,  the  polarity  of  one 
sample  (per  cycle)  will  change  in  (b)  leading  to  a  four  positive,  five 
negative  repetitive  sample  "pattern"  in  (c)  and  a  reference  phase  of  20°. 

More  generally,  when  the  phase  of  the  input  signal  lies  in  the  interval 
[ (n  •  20°  -  10°),  (n  •  20°  +  10°)]  where  n  is  an  integer  or  zero,  the  phase 
of  the  reference  signal  equals  n  •  20°.  Note  that  the  five  positive,  four 
negative  sample  pattern  in  (c)  or  its  negative  -  a  five  negative,  four 
positive  pattern  -  is  the  only  sample  pattern  possible  when  the  input  to 
the  reference  generation  network  is  a  sinusoid  having  a  frequency  fc.  All 
discrete  values  of  the  reference  signal's  phase  are  obtained  from  the 
response  shown  in  (d)  and  its  negative  by  shifting  the  time  origin  in 
increments  equal  to  the  sampling  interval.  The  sample  values  listed  below 
(e)  apply  for  a  unit  amplitude  reference  signal  and  were  obtained  from 
simulation  results. 

When  the  array  output  signal  is  a  sinusoid  and  its  frequency  is 
offset  from  the  center  frequency,  the  phase  of  the  limited  signal  (c)  at 
the  second  filter's  input  increments  in  discrete  20°  jumps  forming  a 
stairstep  pattern  in  time.  The  second  filter  responds  very  rapidly  to 
these  phase  steps  —  on  the  order  of  six  carrier  cycles  --  due  to  this 
filter's  wide  bandwidth.  As  a  result,  the  phase  of  the  reference  signal 
is  essentially  a  stairstep  function  of  time  also. 
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F.  Performance  Measure:  CW  Pattern 


The  performance  measure  selected  for  evaluating  the  simulation  results 
was  the  cw  power  pattern  of  the  array.  The  power  pattern  was  computed  at 
selected  instants  of  time  during  adaption  from  the  set  of  weighting  coeffi¬ 
cients  which  occured  at  those  instants.  The  power  pattern  at  a  frequency 
fc  +  Af  was  defined  as  the  array  output  response  (magnitude  squared  in  dB 
and  phase)  to  a  unit  amplitude  cw  signal  of  frequency  offset  Af  and  variable 
arrival  angle  e: 

(27)  Pz(e,Af)  =  10  log10|Vz(e,Af)|2  exp[j  Arg  VE(e,Af)] 
where  P  (e,Af)  =  cw  power  pattern  phasor 

It 

VyiQibf)  -  array  output  (voltage)  phasor. 

The  expression  for  the  array  output  phasor  in  terms  of  the  set  of 
weighting  coefficients  is  given  in  Appendix  I.  The  computer  simulation 
program  is  listed  in  Appendix  II;  subroutine  REFRN  performs  the  reference 
network  digital  processing. 


SECTION  IV 

COMPUTER  SIMULATION  RESULTS 


A.  Introduction 

The  simulation  parameters  described  in  the  previous  section  were 
assigned  values  as  shown  in  Table  I.  The  parameters  are  listed  in  columns; 
each  line  in  the  table  represents  a  different  simulation  run.  The  parameters 
listed  under  the  heading  "Modifications"  apply  to  modified  computer  simulation 
models  and  will  be  discussed  in  context  as  the  modifications  are  introduced. 
All  parameters  having  dimensions  of  frequency  (or  time)  are  normalized  to 
the  array  center  frequency  fc  (or  its  reciprocal  Tc).  In  all  cases  the 
weighting  coefficients  were  initially  set  to  unity  to  establish  a  broad¬ 
side  pattern.  Because  of  processing  time  restrictions,  it  was  necessary 
to  limit  the  simulation  to  the  case  of  a  two-element  array.  The  important 
response  characteristics  -  especially  those  attributed  to  the  reference  wave¬ 
form  processor  -  do  not  appear  to  be  dependent  on  the  number  of  array 
elements.  These  response  characteristics  are  presented  in  the  following 
sections. 

B.  Effect  of  Reference  Network  Time  Delay 

When  the  code  and  data  rates  are  fixed,  the  selection  of  the  data 
filter  bandwidth  involves  a  tradeoff  between  the  processing  gain  and  the 
time  delay  introduced  by  the  reference  network.  Figure  14  shows  the  effect 
of  time  delay  when  the  data  filter  has  a  bandwidth  equal  to  four  times  the 
data  rate.  The  pattern  magnitude  in  the  desired  signal  direction  (0°)  and 
in  the  interference  direction  (40°)  versus  adaption  time  in  carrier  periods 
(Tc)  are  illustrated.  The  time  delay  results  in  a  180°  phase  difference 
between  the  desired  components  in  the  array  output  and  reference  signals  for 
the  first  fifty  carrier  periods  (50  Tc)  of  each  data  bit  period  (200  Tc) 
since  the  data  is  an  alternating  (worst- case)  sequence .  After  the  initial 
adaptive  transient  reduces  the  desired  signal's  output  level  to  the  correct 
value  (0  dB  pattern  magnitude),  the  character  of  the  signal  response  changes. 
The  feedback  system  attempts  to  reverse  the  sign  of  the  weighting  coeffi¬ 
cients  during  the  delay  error  periods  as  shown  in  Fig.  15.  As  the  weights 
decrease  in  the  process  of  sign  reversal  the  pattern  magnitude  in  the  signal 
direction  decreases.  During  the  remaining  three-quarters  of  each  data  bit 
the  phasing  is  correct  and  the  feedback  system  attempts  to  recover  the  gain 
loss  in  the  signal  direction.  The  pattern  gain  variations  in  the  desired 
signal  direction  can  be  reduced  by  decreasing  the  gain  (bandwidth)  of  the 
feedback  loops  as  shown  in  Figs.  16  and  17.  The  loop  gain  factor  K  has  been 
reduced  by  a  factor  of  five  here.  If  the  loop  gain  is  reduced  even  further 
to  the  point  where  the  response  time  of  the  loops  is  much  larger  than  the 
data  bit  period,  then  the  output  signal  amplitude  is  essentially  constant 
in  the  steady-state.  The  steady-state  amplitude  can  be  calculated  as  shown 
in  Appendix  III  as  the  value  which  minimizes  the  mean-square  error. 
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TABLE  I.  Simulation  Parameter  Assignments 
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The  effect  of  data  modulation  time  delay  in  the  reference  network. 
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reference  network. 


This  steady-state  amplitude  becomes  smaller  as  the  percentage  of  time  when 
bi-phase  errors  are  present  increases,  i.e.  as  the  reference  processor  time 
delay  is  increased.  For  delays  of  one-eighth,  two-eighths,  and  three-eighths 
the  data  bit  period,  the  suppression  of  the  desired  output  signal  below  the 
constant  amplitude  (0  dB)  reference  signal  is  2.5  dB,  6  dB,  and  12  dB 
respectively.  These  figures  apply  for  the  worst-case  (square-wave)  data 
sequence.  For  a  random  data  sequence  the  rate  of  occurence  of  bi -phase 
errors  is  less  by  a  factor  of  two  and  the  corresponding  values  of  sup¬ 
pression  are  1.16  dB,  2.5  dB,  and  4.1  dB,  respectively.  Pattern  adaption 
in  the  interference  direction  in  Figs.  14  and  16  is  very  erratic  as  a 
result  of  the  corrections  caused  by  the  desired  signal. 


It  is  possible  to  equalize  the  reference  processor  delay  by 
employing  wideband  delay  lines  to  delay  1)  the  signal  inputs  to  the  error 
by  signal  multipliers,  2)  the  array  output  signal  applied  to  the  subtractor 
which  forms  the  error  signal,  and  3)  the  local  oscillator  signal  at  the  input 
to  the  second  mixer  in  the  reference  processor  (or  the  corresponding  PN 
sequence  could  be  delayed).  This  approach  allows  other  types  of  data 
modulation  to  be  accommodated,  however,  it  complicates  the  task  of  maintaining 
proper  phasing  conditions  in  the  feedback  loops  of  the  processor.  An 
alternative  to  equalization  which  is  applicable  for  the  case  of  bi -phase  data 
modulation  is  simply  to  inhibit  processing  (i.e.,  hold  the  weights  constant) 
during  periods  when  the  error  signal  is  incorrect.  Figure  18  illustrates 
the  effect  of  holding  the  weights  constant  for  the  first  three-tenths  of 
each  data  bit  period  (T^  =  60  Tc).  This  was  accomplished  in  the  simulation 
program  by  setting  the  error  signal  samples  equal  to  zero  at  the  appropriate 
times:  r 


0 


nf"1  <t  <  nf +  60  Tc 


(28)  E(t)  =  •< 


n  =  0,1,2, ••• 


R(t)-S(t);  elsewhere, 


where  fj  is  the  data  rate.  The  pattern  adapts  smoothly  in  the  desired 
signal  direction  to  the  proper  magnitude  and  the  rejection  of  inter¬ 
ference  is  considerably  improved  compared  to  the  result  shown  in  Fig.  14. 


C.  Effects  of  Code  Synchronization  Error 


When  the  received  signal  code  in  S(t)  and  the  reference  code  are 
not  properly  aligned  in  time,  the  signal  at  the  output  of  the  first  mixer 
in  the  reference  network  is  distorted.  For  small  timing  errors,  the 
smoothed  signal  at  the  data  filter  output  may  be  represented  by  two  com¬ 
ponent  signals:  a  data-carrying  signal  identical  to  that  obtained  when 
no  timing  errors  are  present  except  for  a  reduction  in  amplitude  and  a 
low-level,  noiselike  signal  resulting  from  the  distortion  components  of 
the  codes.  The  bandpass  limiter  removes  the  amplitude  fluctuations  of  the 
composite  signal;  the  phase  fluctuations  due  the  noiselike  component,  how¬ 
ever,  are  retained  at  the  limiter  output  and  therefore  appear  on  the 
reference  signal.  In  the  computer  simulation,  the  phase  of  the  reference 
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signal  was  quantized  so  these  small  phase  variations  were  also  removed.* 
Consequently,  the  reference  signal  in  the  simulation  differed  from  the 
desired  component  of  S(t)  only  in  code  timing  (and  data  modulation  delay, 
of  course). 

The  degradations  caused  by  imperfect  code  timing  and  a  non-zero 
delay  in  the  reference  processing  network  are  expected  to  be  similar 
since  the  reference  signal  is  out  of  phase  by  180°  with  respect  to  the 
desired  reference  signal  over  discrete  time  intervals  in  both  cases. 
However,  the  intervals  are  distributed  differently  in  time  and  their 
durations  are  different  in  the  two  cases.  When  code  timing  is  in  error, 
the  intervals  over  which  the  reference  signal  is  inverted  must  be  much 
shorter  than  the  reference  network  processing  delay  --  i.e.,  on  the  order 
of  one-half  of  a  code  bit  period  or  less  —  if  the  adaptive  array  processor 
is  to  operate  as  intended.  The  rate  at  which  the  inversions  occur  has  an 
average  value  equal  to  one-half  the  code  rate  when  the  code  is  a  (pseudo) 
random  sequence.  The  reference  signal  inversions  caused  by  reference 
network  processing  delay  occur  at  a  much  slower  rate:  one-half  the  data 
rate  on  the  average  when  the  data  bit  stream  is  a  random  sequence.  Code 
timing  errors  equaling  one-quarter  and  one-half  of  a  code  bit  period 
cause  the  responses  to  change  from  those  shown  in  Fig.  18  to  the  responses 
shown  in  Fig.  19.  When  the  error  equals  one-half  the  code  bit  period,  the 
reference  signal  is  inverted  approximately  one-fourth  of  the  time.  This 
is  approximately  the  same  percentage  of  time  the  reference  signal  is 
inverted  as  a  result  of  the  delay  introduced  by  the  reference  processing 
network.  Thus,  it  is  appropriate  to  compare  the  signal  responses  shown 
in  Fig.  14  and  in  Fig.  19  for  the  one-half  code  bit  timing  error  case. 

The  signal  response  in  Fig.  19  is  more  erratic  because  the  reference 
signal  inversions  are  distributed  randomly.  It  is  apparent  from  Fig. 

19  that  the  feedback  circuitry  cannot  change  the  weights  (pattern) 
significantly  over  the  short  inversion  intervals.  Nevertheless  the 
amplitude  of  the  desired  signal  at  the  array  output  is  reduced  relative  to 
the  reference  level  by  nominally  the  same  amount  in  the  two  cases. 

Comparing  the  interference  responses  in  Figs.  18  and  19  shows  that 
the  interference  is  suppressed  less  when  the  code  timing  is  in  error  — 
at  least  for  times  t  less  than  4000  Tc.  In  concept,  the  weights  could 
be  held  constant  from  an  instant  preceding  each  transition  of  the  code 
generated  in  the  array  processor  to  an  instant  following  each  transition 
to  reduce  the  degradation  in  response  caused  by  code  timing  error.  How¬ 
ever,  switching  speed  limitations  of  practical  integrate/hold  circuits 
restricts  the  usefulness  of  such  an  approach  to  cases  where  the  code  rate 
is  relatively  low.  In  many  applications,  the  code  timing  error  should  be 
maintainable  at  less  than  one-fourth  of  a  code  bit  length  and  the  array 
performance  will  be  acceptable  without  resorting  to  holding  the  weights 
constant  during  the  inversion  intervals. 


*This  minor  defect  of  the  simulation  model  is  expected  to  result  in  only 
second-order  effects  since  the  distortion  components  are  usually  small. 
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Fig.  19.  The  effect  of  code  timing  error. 


D.  Effect  of  Desired-Signal  Frequency  Offset 


As  discussed  previously,  phase  alignment  between  the  desired  signals 
in  S(t)  and  R(t)  occurs  when  the  phase  shift  through  the  reference  network 
is  a  multiple  of  360°.  This  condition  is  not  satisfied,  in  general,  when 
the  carrier  frequency  of  the  desired  signal  is  offset  from  the  center 
frequency  since  the  reference  network  phase  shift  is  frequency  dependent.* 

An  analysis  of  the  effects  of  small  frequency  offsets  is  presented  below  and 
is  followed  by  results  of  the  computer  simulation. 


A  few  assumptions  were  made  to  simplify  the  analysis: 


a)  C(t)  =  D(t)  =  1 


b)  An  =  0 


c) 


Af. 


«  1 


(cw  desired  signal) 

(No  interference) 

(small  frequency  offset). 


The  signals  at  the  input  to  the  weighting  coefficients  were  expressed  as 

-  n 

(29)  y2k_6(t)  8  As  sin  uct  +  Aw$t  +  <j>k  -  (1-6)  |- 


k  =  1 ,2  ,•  •  •  ,N 


6  =  0,1 


where  inphase  and  quadrature  signals  from  the  kth  element  correspond  to 
6  equaling  one  and  zero,  respectively.  Signal  phase  angle  <j>.  is  an 
arbitrary  constant.  K 

The  first  step  in  the  analysis  was  to  hypothesize  that  the  array 
output  signal  and  the  reference  signal  can  be  represented  in  the  steady- 
state  as 


(30)  S(t)  =  Ae  sin(w  t  +  Aw  t  +  4>  );  t  -*■  » 
and 

(31)  R(t)  =  Ar  sin[wct  +  Aw^t  +  4> j.  +  y(Au)j.)]  ;  t  -*■  » 


*If  delay  equalization  is  employed  as  discussed  in  Section  IVB,  then  signals 
at  the  error  subtractor  may  be  maintained  in  phase  over  a  range  of  input 
frequencies. 
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where 


Aw  =  array  output  signal  frequency  offset 

<j>  =  array  output  signal  phase 

y(aol>s)  =  reference  network  phase  shift,  (y(0)  =  0) 


AWj,  f  Au^.* 


The  steady-state  error  signal,  given  by  (31)  minus  (30),  is  processed  by 
the  feedback  loops  to  produce  steady-state  weighting  coefficient  variations 
given  by 


(32)  w2k_6(t) 


K  E(t')  y2k-&{V)AV 


- -  (A  sin[(Aw  -Aw  )t  -  <j>_  -  y(Ao)„)  +  v-O-6)  j  ] 

2(A(i)^-Ao)j,)  r£  SEE  E  k  2 


-  Aj.  sin  [(Ao>s  -Ao)j,)t  -  4*2  +  \  “  (1-5)  I"  11 


t  j  k  =  1,  2,  •••,  N  j  6  =  0,1. 

The  weighting  coefficients  modulate  the  input  signals  (29)  resulting  in 
the  steady-state  output  signal 

(33)  S(t)  =  w2k-6^t)  y2k-6(t) 

KA?  N  1 

-  — *— x  ll  {  } 

2  ( Ato  -Aw  )  k=l  6=  0 

S  U 

where 


*Note  that  Awj  and  Aws  are  equal  prior  to  adaption  (t  <  0)  but  are  assumed 
here  to  be  unequal  in  the  steady-state  (t-*»). 
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{  }  =  ^  cos[n>ct  +  Aoi^t  +  4)^  +  y(acos)  ] 

-  jjr  Ar  COS[coct  +  (2  Aa)s  -  Aa>  )t  -  <j>E  -yUo^)  +  2<J>.  -  ( 1  -6 ) 7r  ] 

a 

-  Aj  COS  [oict  +  Ao^t  +  *z] 

+  ^- Aj,  COS  [w^t  +  (2Aojs  -AodjJt  -  <J>2  +  2<j)|^  -  (1  -5 ) it] 

When  this  bracketed  expression  is  summed  over  the  two  values  of  6,  the 

second  and  fourth  terms  cancel  leaving 

1 

(34)  l  {  }  =  A  COS  [w  t  +AO)  t  +  <j>„  +  y(Ao)  )  ] 

6=0  rs,  c  z  2  z 

-  Az  COS  [M  t  +  Au^t  +<j>  ]. 

The  expression  (34)  represents  the  contribution  of  the  kth  element  to 
the  array  output  signal.  It  is  independent  of  the  index  k;  therefore, 
the  element  output  signals  are  all  aligned  in  phase,  i.e.,  they  are 
phase-coherent  in  the  steady-state.  The  cancellation  of  terms  in¬ 
volving  the  input  phase  angles  ^  does  not  occur  if  the  quadrature 
channel  phase  shifts  are  different  from  ninety  degrees.  A  more 
general  analysis  is  necessary  in  this  case  with  anticipated  results 
being  a  degradation  in  signal  coherence  and  the  existence,  in  the  steady- 
state,  of  additional  frequency  components  in  the  array  output  signal. 

In  the  present  case  the  expression  in  (33)  reduces  to  the  form 

NK  A2 

(35)  S(t)  =  - - -  {-A  sin  y(Au  )  sin[w  t  +  Aw_t  +  <J>J 

2(AVAU)E)  ^  CEE 

+  [Ar  cos  y ( Au>e )  -  Ae1  COS  [u>ct  +  Au>Et  +  <J>E]} 

Xt 

by  use  of  a  trigonometric  identity  in  (34).  The  conditions  that  must  be 
imposed  in  order  that  (35)  agree  with  the  original  hypothesis  (30)  are 


(36)  Ae  =  Ar 

& 


COS  y(Au>e) 
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and 


(37) 


-  NK  b  sin  y(au_) 
-i _ Li _ L 

2{Aojs  -AUj.) 


The  first  condition  shows  that  the  amplitude  of  the  output  signal  is 
reduced  below  the  reference  signal  amplitude  by  a  factor  equal  to  the 
cosine  of  the  steady-state  phase  error.  The  two  conditions  may  be  com¬ 
bined  to  obtain  a  transcendental  equation  for  the  output  signal  fre¬ 
quency  offset: 


2 (am.  -  Aui  ) 

(38)  tan  [-  y(a»J]  =  +  -  — 

1  NK  Af 

5  * 

The  two  sides  of  this  equation  are  graphed  in  Fig.  20.  The  straight  line 
representing  the  right  side  of  (38)  has  a  slope  inversely  related  to 
loop  gain  factors  and  shifts  with  constant  slope  as  the  input  frequency 
offset  is  varied.  The  equation  is  satisfied,  i.e.,  the  curves  intersect, 
at  a  value  of  output  frequency  offset  which  is  smaller  than  the  input 
frequency  offset  (aw^Aujs).  The  frequency  offset  of  the  array  output  signal 
cannot  be  reduced  to  zero  since  the  feedback  loops  require  a  steady-state 
error  signal  to  maintain  the  cyclic  variations  of  the  weighting  coefficients. 
If  the  loop  gain  is  small,  the  line  in  Fig.  20  is  nearly  vertical  and 
the  reduction  in  frequency  offset  is  small.  For  increasing  input  fre¬ 
quency  offset,  the  point  of  intersection  moves  up  the  tangent  curve, 
the  steady-state  phase  error  increases  toward  ninety  degrees,  and  the 
output  signal  amplitude  approaches  zero. 

The  computer  simulation  results  given  in  Figs.  21  and  22  confirm 
the  analysis.  The  significant  parameters  obtained  from  these  figures  are 
the  following: 


1)  A  =  0.886  (-1.06  dB) 

u 

2)  Af  -  Af  =  0.00085  f  (weight  oscillation  frequency). 

5  L  C 


These  results  agree  almost  exactly  with  those  computed  analytically  from 
(24),  (36),  and  (38)  using  the  parameter  values  for  this  run,  listed  in 
Table  I.  The  phase  shift  of  the  reference  network  is  28°  as  computed 
from  (36).  Note  that  the  assumptions  of  the  analysis  are  not  completely 
satisfied  here  since  the  desired  signal  contains  bi-phase  PN  coding  and 
low-rate  data  modulation.  The  weighting  coefficients  are  held  constant 
for  a  period  of  time  after  the  data  bit  bi -phase  transition  at  t  =  2000  Tc. 
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Fig.  20.  Results  of  desired  signal  frequency  offset  analysis. 


43 


-90 


3SVHd  NU3JLJLVd 


(8P )  3anJLIN9\M  Nd3J.iVd 


44 


Fig.  21.  The  effect  of  desired  signal  frequency  offset. 
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Fig.  22.  Weighting  coefficient  response  for  result  in  Fig.  21. 


During  this  hold  period,  the  signal  input  and  output  frequencies  are 
identical.  Upon  resuming  processing,  the  cyclic  variation  of  the  weighting 
coefficients  is  restored  as  shown  in  Fig.  22.  When  the  loop  gain  (K)  was 
reduced  by  a  factor  of  four,  the  results  obtained  were  as  expected:  a 
lower  output  amplitude  (-1.95  dB),  a  smaller  weight  oscillation  frequency 
(0.00022  fc),  and  a  larger  phase  error  in  steady-state  (37°).  Figure  23 
shows  simulation  results  for  three  values  of  desired  signal  frequency 
offset:  zero,  0.0005  fc,  and  0.0025  fc.  In  this  result,  cw  interference 
has  been  added,  and  the  data  rate  has  been  increased  significantly.  The 
pattern  response  in  the  desired  signal  direction  is  essentially  unaffected 
by  small  offsets  in  frequency  and  is  reduced  3  to  4  dB  when  the  frequency 
offset  equals  0.0025  fc:  one-half  the  data  rate.  The  addition  of  equal 
amplitude  cw  interference  at  center  frequency  in  this  case  did  not  sig¬ 
nificantly  alter  the  response  to  desired  signal.  The  amount  of  inter¬ 
ference  suppression,  however,  does  appear  in  Fig.  23  to  be  dependent  on 
the  desired-signal  frequency  offset.  The  0.0025  fc  offset  case  in  Fig. 

23  was  repeated  with  data  modulation  removed  from  the  desired  signal 
and  with  continuous  loop  processing  (i.e.,  no  data  "hold"  intervals). 

No  significant  changes  in  performance  were  observed.  The  code  rate  was 
increased  from  0.05  fc  to  0.20  fc  in  a  third  run  and  again  there  was  no 
appreciable  change  in  interference  rejection  over  the  initial  4000  Tc 
adaption  period. 

E  Effect  of  Reference-Network  Processing  Gain 

The  effect  of  waveform  processing  gain  in  the  reference  network 
was  assessed  by  comparing  results  for  different  code  rates  with  the  data 
filter  bandwidth  held  fixed  at  0.02  fc.  Figure  24  illustrates  performance 
for  code  rates  equaling  0.01  fc,  0.05  fc,  and  0.20  fc:  (one -half,  five- 
halves,  and  ten  times  the  data  filter  bandwidth,  respectively).  The 
pattern  response  in  the  interference  direction  has  been  sketched  in  between 
values  computed  at  the  hold  instants  to  make  the  three  cases  more  distin¬ 
guishable.  The  essential  features  to  be  noted  in  Fig.  24  are  the  following 
1)  for  very  small  processing  gain  (the  .01  fc  case),  the  rejection  of  inter 
ference  is  very  erratic,  2)  for  moderate  processing  gains  ( fc ^  =  .05  fc 
and  0.20  fc),  the  response  to  interference  is  much  smoother  with  slightly 
better  rejection  at  the  higher  code  rate.  The  erratic  performance  for 
very  low  code  rates  is  to  be  expected  since  the  data  filter  does  not 
average  the  bi -phase  coded  interference  over  several  bits  of  the  code. 

The  response  to  interference  in  this  case  is  strongly  dependent  on  the 
detailed  structure  of  the  code  sequence  being  used.  Initial  (transient) 
behavior,  in  particular,  is  influenced  by  the  starting  point  in  the  code, 
i.e.  by  the  number  and  distribution  of  transitions  in  the  first  several 
bits  of  the  code.  The  starting  point  for  the  code  employed  in  the  simu¬ 
lation  (see  Fig.  5)  produced  the  rather  surprising  result  of  a  faster 
response  and  more  rejection  with  the  0,01  fc  code  rate  than  for  the  higher 
code  rates  in  the  time  interval  t  <  4000  T~.  However,  it  is  believed  that 
different  code  starting  points  wouTd  result  in  dramatically  different 
initial  responses  for  the  0.01  f  code  rate  case,  some  of  which  would  be 
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Fig.  23.  The  effect  of  desired  signal  frequency  offset 
(with  interference  present). 


much  less  desirable  than  the  response  shown.*  In  addition,  it  is  believed 
that  the  response  for  fed  =  0.01  fc  is  not  uniformly  better  than  the 
responses  shown  for  higher  code  rates  for  t  >  4000  Tc.  Thus,  although 
the  information  presented  is  incomplete  due  to  limits  on  computation  time, 
it  appears  that  the  ratio  of  the  code  rate  to  the  data  filter  bandwidth 
must  exceed  approximately  three  in  a  practical  system. 

F.  Performance  for  Large  Input  Interference 

The  simulation  results  presented  to  this  point  have  been  limited  to 
the  case  where  the  desired  and  interfering  signals  are  equal  in  magnitude. 
The  performance  of  the  processor  for  larger  levels  of  input  interference 
will  now  be  examined.  Before  proceeding  with  the  results,  it  is  desirable 
to  discuss  the  characteristics  of  the  processor  which  affect  its  speed 
of  response.  The  gain  of  the  feedback  loops  (and  thus  their  response  time) 
is  dependent  upon  the  amplitude(s)  of  signal (s)  present  at  the  array  inputs. 
This  dependence  can  be  shown  from  the  defining  equations  (5)  for  the 
weighting  coefficients: 

2N 

(39)  w,  =  K  {yi(t)  (R(t)-  I  'V*)  Vt!l!Baseband 

J  Terms 

Upon  rearrangement,  these  equations  become 
~  2N 

(40)  w.  +  K  {yi(t)J  Baseband  w.  +  K  ^  w.(t)  Cyi(t)  yj(t)}  Ba$eband 

j^i 


=  K  (yi(t)  R(t)}Baseband> 

The  coefficient  multiplying  ith  weight  Wj  in  (40)  is  proportional  to  the 
square  of  (i.e.,  the  power  of)  the  signal  y-j(t)  at  the  input  to  the  ith 
feedback  loop.  By  analogy  with  the  simple,  constant-coefficient  equation 


(41)  w  +  xw  =  y  ;  x,  y  constants, 

whose  solution  is  of  exponential  type,  exp  [-Xt],  it  is  to  be  expected 
that  the  response  time  of  the  ith  feedback  loop  is  strongly  dependent  on 
this  coefficient. 


♦Variation  in  the  response  with  code  starting  point  is  also  expected  at 
higher  code  rates,  although  the  amount  of  variation  should  decrease  as 
the  code  rate  is  increased. 
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The  system  of  equations  (40)  can  be  solved  analytically  for  the 
special  case  of  a  single  cw  signal  incident  upon  the  array  and  a  cw 
reference  signal.  This  is,  of  course,  the  constant  coefficient  case; 
however  the  system  of  equations  is  still  coupled.  Details  of  the  solution 
are  given  in  Appendix  IV.  The  weighting  coefficients  decay  exponentially 
in  this  case  with  a  time  constant  inversely  proportional  to  signal  power: 


This  time  constant  also  appears  in  the  result  which  shows  the  effects 
of  frequency  offset  (38). 

When  interfering  signals  are  present  and  a  processed  reference 
signal  is  used,  the  weighting  coefficient  responses  are  not  necessarily 
given  by  simple  exponentials.  Consequently,  the  notion  of  a  single  time 
constant  as  a  measure  of  convergence  may  not  be  appropriate.  The  above 
results,  however,  indicate  that  the  weighting  coefficients  become  more 
responsive  as  the  power  level  of  either  the  desired  signal  or  interfering 
signals  increases.  This  does  not  imply  that  the  rates  of  change  of  the 
various  signals  in  the  array  output  are  uniformly  affected  by  some  measure 
of  the  total  input  power.  On  the  contrary,  the  different  output  signals, 
or  equivalently,  the  array  pattern  in  different  directions,  may  change 
at  differing  rates.  A  possible  interpretation  here  is  that  the  rate  of 
response  of  the  array  to  a  given  input  signal  is  primarily  determined  by 
the  power  contained  in  that  signal  and  to  a  much  lesser  extent  by  the 
characteristics  of  the  other  input  signals.  That  is,  there  is  some  evidence* 
which  supports  the  conjecture  that  the  array  responds  to  each  signal  more 
or  less  independently  except  when  1)  the  number  of  interfering  sources 
exceed  the  number  of  independent  pattern  nulls  (the  overcons trained-array 
case),  2)  the  angular  separation  between  two  or  more  sources  is  very 
small  or  zero,  and  3)  an  interfering  source  is  sufficiently  large  so  as 
to  produce  non-linear  effects  (saturation  and  limiting)  in  the  array 
circuicry.  Another  situation  leading  to  non-independence  of  desired  and 
interfering  signal  responses  occurs  when  the  response  times  of  the  feed¬ 
back  loops  are  decreased  -  as  a  result  of  increased  interference  power  - 
to  values  less  than  a  code  bit  period.  The  simulation  results  which  follow 
illustrate  this  case. 

Pattern  adaption  results  for  a  20  dB  interference-to-signal  ratio  at 
the  array  inputs  are  shown  in  Fig.  25.  The  element  spacing  was  increased 
to  ten  wavelengths  here  to  illustrate  the  increase  in  angular  resolution 
with  wider  spacing.  The  initial  (broadside)  pattern  is  maximum  (+9  dB)  in 
the  interference  direction  (0°)  and  nearly  zero  (-32  dB)  in  the  desired 
signal  direction  (2.85°).  Corresponding  initial  output  levels  are  +29  dB 


*Results  obtained  by  experimental  testing  with  a  four-element  adaptive 
array  processor  lend  support  to  this  interpretation. 
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Fig.  25.  The  effect  of  processing  loop  bandwidths  approaching  the  code  rate. 


for  interference  and  -32  dB  for  desired  signal:  a  61  dB  interference-to- 
signal  ratio.  The  frequency  offsets  of  both  signals  are  zero.  The  processor 
adapts  quickly  to  reduce  the  large  output  interference  level  to  approximately 
-3  dB  (note  the  -23  dB  pattern  response  in  the  interference  direction. 
Simultaneously,  the  desired  output  signal  increases  from  -32  dB  to  ap¬ 
proximately  -13  dB.  Thereafter,  the  pattern  magnitude  in  the  desired 
signal  direction  increases  slowly  with  a  corresponding  slow  decrease  in 
the  maximum  value  of  the  pattern  in  the  interference  direction.  The  pattern 
variations  in  the  interference  direction  and  the  small  changes  in  the 
weighting  coefficients  in  Fig.  26  represent  the  response  of  the  processor 
to  errors  produced  by  bi-phase  code  transitions  of  the  desired  signal. 

The  essential  characteristics  of  this  response  are  illustrated  by  phasor 
diagrams  of  the  array  output  signals  in  Fig.  27.  The  desired  signal  and 
interference  phasors  are  approximately  equal  in  magnitude  (-6  dB)  at  a 
time  near  t  =  2400  Tc  as  evidenced  by  the  20  dB  pattern  differential  in 
the  two  directions.  Before  a  code  transition,  the  sum  of  output  desired 
signal  and  interference  (nearly)  equals  the  reference  signal  and  the  error 
signal  is  (nearly)  zero.  This  condition  is  represented  by  the  dotted 
phasors  of  Fig.  27a.  Immediately  after  a  code  transition,  an  error  results 
which  is  removed  primarily  by  a  (near)  bi -phase  transition  of  the  output 
interference  phasor  as  shown  in  Fig.  27b.  As  the  output  interference 
phasor  undergoes  this  transition,  its  magnitude  first  decreases  and  then 
increases  as  the  transition  is  completed.  The  "up-and-down"  pattern 
magnitude  variations  in  Fig.  25  in  the  interference  direction  illustrate 
this  behavior  over  many  code  transitions.  Only  a  very  small  change  occurs 
in  the  desired  signal  phasor  of  Fig.  27.  Changes  in  the  array  output 
interference  dominate  over  desired  signal  changes  for  two  reasons: 

1)  the  error  x  input  interference  components  at  the  outputs  of  the 

error  multipliers  in  all  the  feedback  loops  are  much  larger 
(collectively)  than  the  error  x  input  signal  components,  and 

2)  only  small  weight  changes  are  required  to  produce  large  changes 
in  the  output  interference  component  when  the  input  interference 
is  large. 

Thus,  an  error  signal  can  be  eliminated  more  rapidly  by  changing  the  in¬ 
terference  at  the  array  output  rather  than  the  desired  signal.  Full  180° 
phase  transitions  in  the  output  interference  phasor  are  completed  only  when 
several  code  bits  of  the  same  sign  follow  a  code  bit  transition.  Since  this 
code  structure  seldom  occurs  in  pseudorandom  codes,  the  error  signal  is  not 
completely  nulled  between  code  transitions.  Processing  of  the  residual 
error  signal  over  many  code  transitions  results  in  the  cummulative  effects 
observed:  a  steady  increase  in  the  desired  signal  phasor  and  a  cor¬ 
responding  decrease  in  the  interference  phasor  (see  Fig.  27b).  The  rate 
at  which  the  changes  occur  in  Fig.  25  is  strongly  influenced  by  the  ratio 
of  the  feedback-loop  bandwidths  to  the  code  modulation  bandwidth.  The 
behavior  of  the  phasor  diagrams  with  time  in  this  example  indicate  that 
this  ratio  is  approximately  equal  to  one. 
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Fig.  27.  Array  output  conditions:  (a)  immediately  after 
a  code  transition  and  (b)  after  processing. 


A  40  dB  interference-to-desired  signal  input  ratio  produces  the 
pattern  response  shown  in  Fig.  28.  All  parameters  except  interference 
level  are  unchanged  from  the  previous  case.*  For  the  present  case,  the 
phase  of  the  interference  at  the  array  output  is  switched  rapidly  by  180° 
as  a  result  of  changes  in  the  weights  following  a  180°  transition  in  the 
reference  signal's  phase.  Initially,  the  processor  responds  rapidly  to 
reduce  the  output  interference  level  to  approximately  0  dB  (note  the 
-40  dB  pattern  response).  After  the  decay  of  this  initial  transient, 
there  is  no  tendency  for  the  processor  to  increase  the  level  of  the  desired 
signal  at  its  output  to  the  desired  value.  The  loops  respond,  i.e.,  the 
weights  are  changed,  so  that  the  cw  interference  at  the  processor  inputs 
is  bi -phase  modulated  by  the  reference  code.  As  a  result,  the  error  signal 
can  be  made  small  on  the  average  even  though  the  level  of  the  desired  signal 
at  the  output  is  much  smaller  than  the  desired  value.  The  "coded"  output 
interference  in  S(t)  is  processed  by  the  reference  network  to  produce  a 
"clean"  coded  reference  signal.  Since  the  desired  signal  component  at  the 
output  is  down  approximately  30  dB  the  data  modulation  produces  only  small 
amplitude  and  phase  changes  on  the  composite  output  signal  S(t),  i.e.,  the 


*This  level  of  interference  (An=100)  is  close  to  the  maximum  permissible 
value  for  the  chosen  gain  constant  and  sampling  rate.  Larger  levels  of 
interference  produce  instability  in  the  digital  feedback  loops  of  the 
simulation  (see  Sec.  Ill .0) . 
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The  effect  of  processing  loop  bandwidths 
exceeding  the  code  rate. 


data  is  essentially  removed  in  this  case.  To  determine  if  this  behavior 
is  unique  to  the  case  where  the  frequency  of  the  cw  interference  equals  the 
center  frequency  of  the  processor,  the  frequency  of  the  interference  was 
offset  by  0.01  fc  and  the  simulation  repeated.*  The  pattern  magnitude 
responses  remained  essentially  unchanged;  however,  the  pattern  phase  angle 
in  the  interference  direction  decreased  linearly  except  at  code  transition 
instants  where  a  180°  step  change  was  observed.  This  result  implies  that 
the  processor  tends  to  remove  the  frequency  offset  of  the  interfering  signal 
in  addition  to  bi -phase  modulating  it  with  the  code  rather  than  responding 
properly  to  the  much  smaller  desired  signal. 

The  pattern  response  when  the  initial  pattern  is  broadside,  the 
desired  signal  arrives  at  an  angle  of  0°,  and  the  interference  at  40°  is 
shown  in  Fig.  29.  The  element  spacing  here  is  one-half  wavelength.  The 
pattern  response  in  the  direction  of  the  desired  signal  reduces  from  9  dB  to 
6  dB  as  a  40  dB  null  is  formed  in  the  interference  direction.  At  the  array 
output,  the  desired  and  interfering  signals  have  steady-state  magnitudes  of 
two  and  one,  respectively.  They  are  180°  out-of-phase  except  during  data 
hold  intervals  when  processing  is  interrupted.  Bi -phase  coding  of  the 
output  interference  component  occurs  as  in  the  previous  example.  In  this 
case,  however,  the  bi-phase  data  modulation  is  retained  on  the  composite 
array  output  signal  since  the  desired  signal  component  at  the  output  is 
larger  than  the  interference. 

It  is  apparent  from  the  results  shown  in  Figs.  28  and  29  that  the 
bandwidths  of  the  processing  loops  have  increased  to  the  point  where  they 
greatly  exceed  the  code  modulation  bandwidth  as  a  result  of  increasing  the 
power  in  the  interfering  signal.  This  condition  occurs  in  the  simulation 
because  the  analytically-derived  feedback  loops  of  Fig.  2  contain  no  band- 
limited  elements**  to  limit  loop  response  time  as  the  levels  of  the  input 
signals  are  increased.  An  unlimited  bandwidth  in  the  ideal  processor  is, 
of  course,  consistent  with  the  theoretical  objective:  to  minimize  the 
squared  error.  As  the  processor  becomes  more  responsive  with  higher 
signal  levels,  the  composite  array  output  signal  is  changed  more  quickly 
to  make  it  conform  to  the  reference  signal  and  smaller  errors  occur  as  a 
result.  The  responses  shown  in  Figs.  28  and  29  confirm  this  behavior: 
the  error  signal  is  maintained  at  a  very  small  value  throughout  the 
processing  interval.  It  should  be  noted  from  these  figures,  however,  that 
the  output  interference  is  not  reduced  to  a  small  value;  also,  the  amplitude 
error  of  the  desired  signal  is  not  reduced  to  zero.  The  important  conclusion 


*The  phase  shift  through  the  reference  network  at  this  frequency  (1.01  fc) 
is  approximately  200°.  The  phase  shift  equals  0°  when  the  frequency  offset 
equals  zero. 

**The  ith  error  multiplier  output  is  the  product  E(t)  y-j (t)  rather  than  the 
convolution  of  this  product  with  a  low-pass  impulse  response  as  would  be 
the  case  for  an  error  multiplier  having  a  limited  output  bandwidth. 
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effect  of  processing  loop  bandwidths  exceeding  the  code  rate. 


to  be  drawn  from  these  results  is  that  minimization  of  the  error  signal  is 
equivalent  to  minimization  of  interference  and  proper  adjustment  of  desired 
signal  at  the  array  output  only  when  the  loop  bandwidths  are  less  than  the 
code  modulation  bandwidth.  This  conclusion  is  consistent  with  the  discussion 
in  Sec.  IIB  regarding  the  correlation  properties  of  the  input  and  reference 
signals  required  for  the  array  to  operate  properly.  It  was  noted  there  that 
the  spectral  width  of  the  reference  signal  should  greatly  exceed  loop  band- 
widths  so  that  the  correlation  between  the  reference  signal  and  the  inter¬ 
ference  over  time  intervals  having  a  duration  equal  to  the  reciprocal  of 
the  loop  bandwidth  is  small. 

The  fact  that  loop  bandwidths  in  the  simulation  become  excessively 
large  as  the  interference  level  is  increased  does  not  necessarily  imply 
that  such  will  be  the  case  in  practical  array  processors.  A  fairly  large 
value  was  used  for  loop  gain  constant  K  in  the  simulation  to  restrict 
program  run  time;  this  value  may  not  be  achievable  in  a  practical  processor. 
Moreover,  a  non-zero  error  multiplier  rise  time  and  integrator  slew-rate 
limitations  may  well  impose  upper  limits  on  processing  (loop)  bandwidths 
which  are  less  than  the  code  modulation  bandwidth.  The  simulation  results 
for  high-level  interference  would  not  apply  in  this  case  since  the  simu¬ 
lation  model  is  incomplete.  The  limitations  in  dynamic  range  of  practical 
array  processors  must  also  be  considered  if  an  extreme  interference  con¬ 
dition  is  to  be  tolerated.  The  occurrence  of  amplifier  saturation,  inter¬ 
modulation  product  generation,  and  non-linearities  in  the  feedback  circuits 
will  undoubtedly  alter  the  responses  from  those  given  in  the  simulation. 

The  simulation  results,  however,  do  illustrate  the  performance  characteristics 
expected  from  an  array  processor  ideally  implemented  in  the  minimum  mean- 
square  error  configuration.  In  the  following  section,  the  possibility  of 
modifying  the  processor's  configuration  to  restrict  the  loop  bandwidth  in 
the  simulation  for  the  high-level  interference  case  is  examined. 

G.  Results  for  Modified  Simulation  Models 


1.  Loop  Gain  Control  (AGC) 

A  time-varying  loop  gain  factor  G(t)  was  introduced  into  each 
feedback  loop  in  the  simulation  to  determine  the  feasibility  of  automatic 
gain  control  (AGC).  Factor  G(t)  had  an  initial  (t=0)  value  of  unity  and 
decayed  exponentially  to  a  final  value  equal  to  the  ratio  of  desired 
signal  power  to  total  input  power;  A£/(A|+Afi).  Figure  30  shows  the  results 
for  40  dB  input  interference  to  signal  ratio  and  a  G(t)  decay  time  constant 
of  250  Tc.  The  response  for  t  <  400  Tc  is,  as  expected,  very  similar  to 
the  response  shown  in  Fig.  29.  As  the  loop  bandwidths  become  smaller  than 
the  code  bandwidth  due  to  the  decay  of  G(t),  the  interference  is  increasingly 
rejected  at  the  array  output.  However,  the  pattern  magnitude  in  the  de¬ 
sired  signal  direction  does  not  change  significantly  indicating  that  a  40  dB 
lnop  gain  reduction  to  achieve  satisfactory  interference  rejection  results 
in  an  unsatisfactory  (slow)  response  to  desired  signal.  This  result  is  to 
be  expected  since  both  desired  and  interference  components  at  the  error 
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Fig.  30.  The  effect  of  loop  gain  reduction  under  high-level  interference  conditions. 


multiplier  outputs  are  reduced  equally  by  G(t).  Although  the  gain  control 
is  in  the  feedback  paths  here,  essentially  equivalent  results  are  anticipated 
with  input  (array  element)  AGC  methods  since  the  gain  reduction  is  again 
common  to  both  desired  signal  and  interference.  The  results  in  Fig.  30  sug¬ 
gest  that  AGC  methods  might  be  effective  provided  1)  the  initial  level  of 
desired  signal  in  the  array  output  is  adequate,  i.e.,  the  pattern  does  not 
have  a  null  in  the  direction  of  desired  signal  initially,  and  2)  the  rapid 
changes  in  the  pattern  occuring  initially  when  the  interference  null  is 
formed  do  not  produce  a  null  in  the  direction  of  the  desired  signal.  The 
use  of  a  moderate  amount  of  gain  control  in  the  array  element  input  cir¬ 
cuitry  might  be  desirable  to  prevent  non-linearities  (saturation)  in  the 
weighting  coefficient  multipliers  when  the  input  signal  levels  are  large. 

2.  Error-Multiplier  Output  Non-linearity 

The  effect  of  amplitude-limiting  (saturation)  in  the  error  multi nliers 
under  high-level  interference  conditions  was  investigated  in  the  simulation. 
The  output  samples  E ( t j )y^* (t j )  of  the  error  multipliers  were  adjusted  in 
amplitude  in  accordance  with  the  non-linear  characteristic  shown  in  Fig.  31. 


Fig.  31.  Simulation  characteristic  for  error  multiplier 
output  amplitude  saturation. 
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For  small  inputs,  the  characteristic  is  linear  and  the  gain  equals  one. 

For  large  inputs,  the  output  amplitudes  are  compressed  according  to  a 
logarithmic  function.  Simulation  results  for  a  40  dB  input  interference  to 
signal  ratio  are  shown  in  Fig.  32.  These  results  show  that  the  non¬ 
linearity  effectively  reduces  loop  gain  initially  when  the  output  inter¬ 
ference  (and  thus  the  error  signal)  is  very  large.  The  rates  of  pattern 
response  shown  for  t  <  200  Tc  are  considerably  smaller  than  those  without 
the  non-linearity  (see  Fig.  29).  As  the  output  error  becomes  small,  how¬ 
ever,  and  linear  operating  conditions  are  restored,  the  gains  (and  bandwidths) 
of  the  feedback  loops  increase.  For  t  >  200  Tc,  the  character  of  the 
response  is  identical  to  that  obtained  previously,  i.e.,  bi -phase  switching 
of  the  interference  occurs  to  remove  errors  resulting  from  the  reference- 
code  phase  transitions.  The  responsiveness  of  the  feedback  loops  when  the 
error  is  small  is  not  affected  by  the  non-linearity. 
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SECTION  VI 
SUMMARY 


A  receiving  array  capable  of  adaptive  pattern  shaping  has  been 
described.  The  signals  at  the  outputs  of  the  array  elements  are 
(effectively)  adjusted  in  amplitude  and  phase  by  a  feedback  control 
system  designed  to  minimize  the  difference  between  the  array  output  signal 
and  a  reference  signal.  The  reference  signal  is  derived  from  the  array  out¬ 
put  signal  by  a  waveform  processing  network.  The  particular  waveform  processor 
employed  accepts  digital  communication  signals  which  are  bi-phase  modulated 
with  the  modulo-2  sum  of  a  code  and  a  random  data  sequence. 

The  initial  (transient)  adaption  performance  of  a  two-element  adaptive 
array  processor  has  been  simulated  on  a  digital  computer  for  the  case  of  cw 
input  interference.  The  simulation  results  show  that  the  waveform  processor's 
time  delay  (or  equivalently,  its  bandwidth)  is  a  significant  parameter  af¬ 
fecting  performance.  The  waveform  processor  delays  the  desired  signal's  data 
modulation  and,  in  addition,  shifts  its  carrier  phase  by  an  amount  dependent 
on  the  carrier  frequency.  These  alterations  reduce  the  correlation  between 
the  reference  signal  and  the  desired  signal  at  the  array  output.  As  a  con¬ 
sequence,  the  desired  signal  is  reduced  in  amplitude  (i.e.,  partially  sup¬ 
pressed)  at  the  array  output  to  a  level  below  that  of  the  constant-amplitude 
reference  signal.  The  worst-case  suppression  due  to  data  modulation  delay 
is  shown  to  be  approximately  2.5  dB,  6  dB,  and  12  dB  for  time  delays  of  one- 
eighth,  two-eighths,  and  three-eighths  of  the  data  bit  period,  respectively. 

The  suppression  factor  due  to  waveform  processor  phase  shift  (with  no  data 
modulation  present  on  the  desired  signal)  is  equal  to  the  cosine  of  the  phase 
shift.  The  simulation  results  show  that  the  weighting  coefficients  undergo 
periodic  variations  in  the  steady-state  in  an  attempt  to  shift  the  carrier 
frequency  of  the  desired  signal  to  *  “nlue  where  the  waveform  processor's 
phase  shift  is  a  multiple  of  36^,  <ne  phase  error  is  reduced  by  an  amount 
dependent  on  the  gain  of  the  fueauack  loops.  Complete  elimination  of  phase 
error  is  not  possible  sire-,  a  residual  error  signal  is  necessary  to  force 
the  weighting  coefficients  to  exhibit  e  periodic  behavior. 

The  delay  introduced  by  the  waveform  processor  also  influences  the 
interference  reject un  capabilities  of  the  array.  The  simulation  results 
show  that  suppression  o-p  desired  signal  during  the  initial  adaption  period 
is  accompanied  by  a  decrease  in  interference  rejection.  A  technique  which 
compensates  for  data  delay  errors  by  processing  only  during  appropriate 
time  intervals  was  described  and  its  effectiveness  demonstrated.  A  delay 
equalization  technique  which  could  be  used  to  improve  the  processor's 
performance  was  also  described. 
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The  correlation  between  the  reference  signal  and  desired  signal  is 
also  reduced  by  imperfect  synchronization  of  the  code  generated  within  the 
processor.  The  simulation  results  show  that  timing  errors  of  one-quarter 
and  one-half  a  code  bit  period  result  in  approximately  3  dB  and  6  dB  signal 
suppression,  respectively.  The  results  indicate  that  the  "steady-state" 
performance  of  the  array  processor  is  relatively  independent  of  the  wave¬ 
form  processor's  processing  gain  -  provided  reasonable  processing  gains 
(code  rates)  are  employed.  Increasing  the  code  rate  does  improve  the 
initial  rate  of  interference  rejection,  however. 

The  simulation  results  show  that  the  bandwidths  of  the  feedback 
loops  in  the  array  processor  increase  as  the  power  in  the  input  signals 
is  increased.  Unsatisfactory  performance  was  observed  when  the  input 
interference  power  levels  were  increased  to  the  point  where  processing 
loop  bandwidths  exceeded  the  bandwidth  of  the  pseudonoise  code  modu¬ 
lation.  The  array  processor  responds  only  to  the  interfering  signal 
in  this  case  and  alters  its  amplitude  and  phase  modulation  to  match 
the  reference  signal.  The  level  of  the  desired  signal  at  the  array 
output  under  these  conditions  is  shown  to  depend  on  the  initial  values  of 
the  weighting  coefficients.  The  performance  of  practical  array  processors 
under  high-level  interference  conditions  is  expected  to  differ  from  the 
simulation  results  because  of  limitations  in  the  bandwidth  and  dynamic 
range  of  components  used  in  the  implementation. 

The  introduction  of  automatic  gain  control  into  the  feedback  loops 
of  the  simulation  model  resulted  in  smaller  processing  bandwidths  and 
satisfactory  interference  rejection  but  unsatisfactory  (slow)  response  to 
desired  signal.  It  was  shown  that  amplitude  limiting  in  the  error  multi¬ 
pliers  under  large  (initial)  error  signal  conditions  reduces  the  processing 
bandwidths  only  until  the  error  becomes  small. 

The  simulation  results  presented  in  this  report  provide  basic  in¬ 
formation  regarding  the  transient  performance  of  adaptive  array  processors 
having  a  waveform-processed  reference  signal.  The  coverage,  however,  is 
incomplete  in  several  respects.  Further  investigation  is  required  to  de¬ 
termine  how  processor  performance  degrades  as  the  angular  separation  be¬ 
tween  the  sources  of  the  desired  and  interfering  signals  is  reduced  to 
small  values,  i.e.,  to  determine  the  spatial  resolution  properties  of  the 
array.  Processor  performance  when  receiver  noise  is  non-negligible  and 
when  the  interfering  signal  has  a  non-zero  spectral  width  should  also  be 
determined.  Currently,  experiments  are  being  conducted  to  obtain  in¬ 
formation  to  augment  the  simulation  results  contained  in  this  report. 

An  experimental  four-element  adaptive  array  processor  has  been  instrumented 
and  initial  feasibility  tests  performed.  The  data  to  be  obtained  from  the 
experiments  will  provide  a  more  complete  basis  for  designing  practical 
adaptive  array  processors. 
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APPENDIX  I 

CW  PATTERN  COMPUTATION 


The  array  geometry  and  notation  of  Fig.  4  were  used  in  computing 
the  cw  pattern.  Antenna  elements  having  ideal  (isotropic)  patterns  were 
assumed;  the  received  signal  from  each  element  was  assumed  to  be  a  unit 
amplitude  sinusoid  of  frequency  f  +  Af.  The  array  output  signal  was 
represented  as 

.  ,  _  j(w  +Aw)t 

(46)  Vs(t,e,Af)  =  Re{  Vs(e,Af)  e  c  } 

where  the  array  output  (sum)  phasor  is  expressed  as 

> 

)  . 


(47)  y<Mf)  .  (tfn  eM"<9>  +  W.n 


Each  term  in  the  summation  corresponds  to  the  phasor  contribution  of  a  pair 
of  elements  spaced  equidistant  from  the  array  phase  center.  Coefficient  ~ 
tfn  is  the  complex  weighting  coefficient  at  frequency  fc  +  Af  of  the  nth  element 
to  the  right  of  the  array  phase  center;  phase  <j>n(e)  represents  the  phase  shift 
due  to  path  delay  to  that  element.  When  (47)  is  expressed  in  the  notation 
used  in  Sec.  II IB  the  result  is  as  follows: 


(48)  Vz(e,Af) 


I  T  w  ,  -  w  sin  w  cosf  \ 

tf-2,6,10,-.  [L  £  V2fc/  £  V2fcyJ 

dXc  0  +  1 


T  +1 


sin  e 


w 


£+1 


'  V2Sin(lf)-  i  V2cos(tf) 

’Hc(1  + 


,-n¥  +i 


sin  by. 
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APPENDIX  II 

COMPUTER  SIMULATION  PROGRAM 


FTM»  9j  L»  T  - - 

C  CODED  ADAPTIVE  ARRAY  SIMULATION 

C 

C 

C  SPC=NO.  OF  SAMPLES  PER  CYCLE  OF  CARRIER  AT  CENTER  FREQ. 


C 

C 

C***ARRAY 

C 


PROCESSOR  PARAMETERS? 


MUM=MO.  OF  ELEMENTS  ( EVEN  INTEGER) 

DI  S=  ELEMENT  SEPARATION  IN  WAVELENGTH S  AT  CENTER  FREQ. 
GAIN=RATIO  OF  LOOP  GAIN  (K)  TO  ARRAY  CENTER  FREQ • 

W(  I  )  =  CHANNEL  I  WEIGHTING  COEFFICIENT 


C 
C 
C 
C 
C 

C*** INPUT  SIGNAL  PARAMETERS; 


C 

C  THS= SIGNAL  ARRIVAL  ANGLE  FROM  BROADSIDE  (SPATIAL  DEG.) 

C  THM= INTERFERENCE  ARRIVAL  ANGLE 

C  A  S= S I GN  AL  AM PL I TUDE 

C  AM=I NTERFERENCE  AMPLITUDE 

C  FFS=  SI GNAL  FRACTIONAL  FREQUENCY  OFFSET 

C  FFM=I NTERFERENCE  FRACTIONAL  FREQUENCY  OFFSET 

C  ALS0=SI GNAL  INITIAL  PHASE  AT  ARRAY  PHASE  CENTER  (DEG.) 

C  ALMO= INTERFERENCE  INITIAL  PHASE  AT  ARRAY  PHASE  CENTER 

C  CODR=RATI 0  OF  CODE  CLOCK  FREQ.  TO  ARRAY  CENTER  FREQ. 

C  DATR=RATI 0  OF  DATA  CLOCK  FREQ.  TO  ARRAY  CENTER  FREQ. 

C 


C*« PREFERENCE  SIGNAL/NETWORK  PARAMETERS; 

C 

C  IDREF=TYPE  OF  REFERENCE  SIGNAL  GENERATION  DESIRED 

C  1.  FIXED 

C  8.  LINEAR  WAVEFORM  PROCESSING 

C  3.  LIMITED  WAVEFORM  PROCESSING 

C 

C  BWF1=RA  TI 0  OF  5.3-DB  BANDWIDTH  OF  FILTER  1  TO  CEN.  FREQ. 

C  OWF2=  RATI  0  OF  5.3-DB  BANDWIDTH  OF  FILTER  2  TO  CEN.  FREQ. 

C  CDOFF=RATI 0  OF  LOCAL  CODE  TIME  DELAY  TO  CODE  CLK.  PERIOD 

C  HOLD= NON- PROCESSING  INTERVAL  DURATION  IN  CARRIER  CYCLES 

C  AR= REFERENCE  AMPLITUDE  ( I DREF= 1 ) 

C  FFR= REFERENCE  FRACTIONAL  FREQ.  OFFSET 

C  ALRfl  =  REFERENCE  INITIAL  PHASE  (DEG.) 

C  GREF= TOTAL  REFERENCE  NETWORK  GAIN  ( I DREF=  2) 

C  ARL=LIMI  TF.D  REFERENCE  AMPLITUDE  ( I  DREF“3) 

C 

C*»«- pro  GRAM  CONTROL  PARAMETERS; 

C 

C  JWGHT=NO.  OF  SAMPLES  PER  WEIGHTING  COEFF.  PRINTOUT 

C  .JCNT=NO.  OF  SAMPLES  PER  STOP/CONTINIJE  CONTROLS 

C  JPTRN=NO.  OF  SAMPLES  PER  PATTERN  PRINTOUT 

C  PFF=FRAC  TI ONAL  FREQ.  OFFSET  FOR  PATTERN  COMPUTATION 

C  PTNI N=PATTERN  MINIMUM  ANGLE 

C  PTMAX=PATTERN  MAXIMUM  ANGLE 

C  PTI NC=PATTERN  INCREMENT  ANGLE 

C 
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c 

c 

c 

PROGRAM  CAAS 

DIMENSION  WC 16)*YC 16)*  FRKC i 6) * PHSGC 8 ) * PHMC 8 ) > I CC 10) 
COMMOM  FJ*  I  DREF*  JHLD*  JBLAY*  f  ^>C*PI*  CODR 
COMMON  AR* ARC*  CODLO* SUM*  REF* ERR 
900  WRI TEC 1*10) 

10  FORMATC  24H  EMTER  SPG  MUM  DI S  GAIN) 

REA  IK  1,  *  )SPC*  MIJM,  DI  STRAIN 

MUM l=MUM- 1 

MUM 2= 2* MUM 

TK=  GAI M/ SPC 

WRI TEC 1* 1  A) 

1  A  FORMATC  23H  EMTER  INITIAL  WEIGHTS) 

READC  1  *.*  >  C  WC  I  >  »  I  =  1  >  M  UM  2  > 

WRI TEC 1* 16) 

16  FORMATC 21H  EMTER  AS  AM  THS  THM) 

READC  1  *  *  )  AS*  AM*  THS*  THM 
WRI TEC 1* 18) 

IB  FORMATC 25H  EMTER  FFS  FFM  ALSO  ALMO) 

READC 1  *  *  )  FFS*  FFM*  AL SO j  ALMO 
PI=3. 1415926 

PD1 =PI *  C  1  .  +  FFS)*DIS*SIMCTHS»PI/180.) 

PD2=PI*C 1 •  +  FFM ) *DI S*  SI NC  THM* PI /ISO*) 

DO  20  K= 1 »  MUM  1  *  2 
ZK=K 

PHSGC  X)  =  £><*PD1  +  ALSOaPI  / 180* 

PHSGC  K+1)  =  -?:K*PD1  +  ALS0*PI/ 180* 

PHMC  X ) =£KaPD2+ALM0*PI / 180 • 

20  PHMCK+1 )=-ZH»PD2+ALM0*PI/180* 

Ql=PIaC l*+FFS)/2* 

Q2=PI *  C  1  *  +  FFM)/2* 

PSI=?.aPI*C 1 *+FFS)/SPC 
PM  I  =  2 • aPI *  C 1 .+ FFM) /SPC 
WRITEC 1*  22) 

22  FORMATC 17H  EMTER  CODR  DATR) 

READC  1  *  *  )  CODR*  DATR 
JCMAX=  SPC/C9DR+ • 5 
JDM4X- SPC/DATR+ • 5 
JCQDE= JC'MAX 
JDATA= JDMAX 
DO  2 A  1=1, 4 
ICC  I )  =  1 

2 A  ICC  I  +  4) =0 

I DATA= 1  • 

I DT=- 1 
I  CODL=  1 

70  WRITEC  1*72) 

72  FORMATC 24H  EMTER  JWGHT  JCMT  JPTRM) 
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o  o  o 


PEADC  1,*)  JWGHT,  JCNT,  JPTRN 
WRI  TE  C  1 , 7  4  ) 

74  FORMATC  29H  EN3 TER  PFF  PTMIM  PTMAX  PTIJJC) 

READC  I,  *  >  PFF,  PTMIN,  PTM4X,  PTIMC 
PPQ=PI DI S*C  1 .  +  PFF) 

ESN=SIMCPI  *PFF/2»  ) 

ECS=C0SCPI*PFF/2. ) 

KK-M=(  CPTMAX-PTMIN)/PTINC)+1 . 5 
JW=0 
JC>3=0 
JPAT=0 
FJ=- 1 . 

CALL  REFRN 
CALL  PHSC  2,  1  2) 

WRI  TE(  1  ,  80)  SPC,  MUM*  DI  S,  GAI  M,  THS..  THM,  AS,  AM,  FFS>  FFM,  ALSO, 

CALMO, CODR, JCMAX, DATR, JDMAX,  PFF 

80  FORMATC  4HSPC=,  F5«  I , /,  4HMIJM=,  1 2,  10X,  4HDI  S=>  F7.  3,  5X,  5HGAIM=, 
CE10.4, /, 4HTHS=, F7.3, 5X, 4KTHM=, F7.3, /, 3HAS=,F8.4, 5X, 3HAM=, 

CF8 . 4,  /,  4HFFS=,  E 1 0 . 4,  2X,  4HFFM=,  El  0 . 4,  /, 5HAL S0=, F7. 3, 4X, 
C5HALM0=, F7. 3, /, 5HC0DR=, El 0. 4, 1 7X, 6HJCMAX=, I  6, /, 5HDATR=, 

CE10.4,  1  7X,  6HJDMAX=,  I  6, /,  4HPFF=,  E10.4, ///,  “INITIAL  WEIGHTS?  _") 
CALL  PHSC2,NtJM2) 

WRI TEC 1,90>CWCI),I=1, NUM2) 

90  FORMATC  32C  El  2*  5,  /,  1  7X)  ) 

118  CALL  PHSC  4,0) 

CALL  PHSC  2,  3) 

WRITEC 1, 100) 

100  FORMATC 16X, 6HWC J+ 1 ) , 8X, 6HSUMC J) , 9X, 6HREFC J) , 9X,  6HERRC J),//) 

ihhhhhhhhhhkjTART  Qp  Loop 

120  IFC JCODE- JCMAX) 160, 125 
125  JCO  DE=0 
JDL=0 

IC0DE=ICC6)  +  ICC7) 

DO  130  1  =  2,7 
1 1  =9- 1 

130  ICCII)=ICCII-1) 

IFC I  CODE- 1)  140,  140,  135 
135  I CQDE=0 
140  ICC  1 )  =  I  CODE 

IFC JDATA-JDMAX)  1  50,  145 
145  JDATA=0 
JHL  D=  0 

I DATA= I DATA+ I DT 
I DT=- I DT 

150  I  MDD=  I  CODE+  I  DATA 
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IFC  I  MOD"  1)1  55*  1  57#  155 

1 55  ASC=AS 

GO  TO  (  1 56# 160# 1 60># I DREF 

156  ARC=AR 

GO  TO  190 

157  ASC=-4S 

GO  TO  C 1 58# 1 60# 1 60)#  1 DREF 

158  ARC=-AR 
GO  TO  190 

160  IFC  JDLAY-JDL)  170#  165#  1  70 

165  I  CO  DL  =  I  CO  DE 
170  IFCICODL)  175#  175,  180 
175  CODLO"  1  • 

GO  TO  190 
1 80  C0DL0=-1. 

190  DO  000  K=  1  #  MIJM 
I  =  2»K 

YC I  - 1  ) = ASC* SIM ( PHSGC  X)  ) + AM*  SI MC  PHMCK)  ) 

YC I )=ASCttSIMCPHSGCK)”01  >  +  AM»SIiMC  PHMCK) -02) 

PHSGCK)=PHSGCK)+PSI 
I FC  PHSGC  X ) “ 6«  283 1 853)  19  5#  192 
192  PHSGCK)=PHSGCK>“6. 2831853 
195  PHMCK) =PHMC  K)+PMI 

IFCPHMCK)-6. 2831853)200# 193 
198  PHMCK)=PHMCK)-6»2831853 

200  COMTIvMUE 
SUM=0. 

DO  20  5  I  =  1  #  M  UM  2 
20  5  SUM” SUM+Y  C I )*WC I ) 

CALL  REFRM 

320  DO  330  1  =  1#  MUM 2 
FBKCI )=YCI )*ERR 
330  WCI)  =  WCI  )  +  T.K*F8KC  I  ) 

IFC  -JW- JWGHT)  350#  340 
340  JW=0 

CALL  PHSC2#  1) 

WRI  TEC  1  #  342)  FJ#  WC 1 ) #  SUM#  REF#  ERR 
342  FORMATC  2H.J=#  F6«  0#  4C  3X#E12.5)) 

DO  344  1  =  2#  MUM 2 
CALL  PH SC  2# 1 ) 

WRI TEC  1#  346)WCI ) 

344  CONTI MUE 

346  FORMATC 1 IX# El  2. 5) 

350  IFC JPAT-UPTRM) 370# 500 
370  IFC JCM-JCMT)4Q0# 380 
380  CALL  PHSC  2#  1  ) 

WRI TEC  1 #  384) 

384  FORMATC "CONTINUE’  EMTER  1=YES  0=NO  _  > 
READC  1#*)  JDEC 
IFC.JDF.C)  799#  799#  390 
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390  JCN=0 

WRI TEC 1>  72) 

READC  l,  *) JWGHT*  JCMT* -JPTRM 
-400  F«J=FJ+1. 

JC0DE=JC3DE+ 1 
<JDATA= JDATA+  1 
JDL= JDL  +  1 
JHLD= -JHLD+  1 
JW=-JW+1 
JCM=JCM+1 
<JPAT= JPAT+  i 
GO  TO  1 20 
C 

C»»*»*****»»ft»*»*»**E'JD  OF  LOO?  HERE?  PATTERN  PROG.  BELOW***** 
C 

500  JPAT=0 

PTH=PTMIM 
CALL  PH. SC  3>  2) 

CALL  PHSC  2>  1 ) 

WRI  TEC  1 >  502)  FJ 

502  FORMATC  2H.J=#  F6. 0.  3X,  7HPATTERM) 

DO  530  K*<=  1  >  KXM 

SM  R=  0  • 

SM  I  =  0  • 

»UM4sMUM8-8 

DO  510  K=  P.i  MtJM4>  A 
I  =  C  X- 2) /2+ 1 
DO  510  L=l>  3,2 
LL=2~L 
X7,K=I*LL 

7 1=^X*PP0*SIUCPTH»PI/150.  ) 

7  1  R=COSC?.  1  ) 

?:1Q=SIMC7.1) 

LK=K-LL 

RR=WCLK)-WCLK+1 )*ESN 
QQ=WCLK+1 )*ECS 
SMR=SMR+X  lR*RR+7, 10*00 
510  SMI  =  SMI +7.  lf)»RR-7.  1R*Q0 

PO  WR=  4 . 34 29  5*  AL 0  6 C  SM R»  SMR+  SM I  * SM I ) 

ARGU=  SMI /SMR 

AMGH=  57  »  29  58*  ATAM  C  ARGIJ) 

CALL  PHSC  2#  1 ) 

WRI TEC  1 »  525) PTH j  POUR*  AMGH 
525  FORMATC  3X.»  F7»  2/  5X.»  F9  •  4.»  5X*F9.4) 

530  PTH=PTH+PTI  NIC 
CALL  PHSC  3>  2) 

GO  TO  370 

799  CALL  PHSC  5*0) 

WRI TEC  1*800) 

800  FORMATC 32H  AMOTHER  RDM  EMTER  1=YES  0=M0) 

READC 1,*).JG0 

IFC JG0)999, 999/900 
999  STOP 
EMD 


C 
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o  a  o  a  o  o 


USE  Ivl  TH  PROGRAM  CAAS:  CODED  ADAPTIVE  ARRAY  SIMULATION 


25 

26 


30 

32 


AO 


A2 


AA 


1 


SUBROUTINE  REFRN 

DIMEM  SION  0 1 C  A)  *  0 1 T!  A)  *  0 1 TTC  A)  *  Bl !  A)  #  EXCO  1  (  A; 
DIMENSION  02! A) #  02T! A) #  02TT!  A)  #  B2!  A)  #  EXC02!  A) 
COMMON  FJ#  I  DREF#  JHLD#  -JDLAY#  SPC#  PI  #  CO  DR 
COMMON  AR#  ARC* CODLO#  SUM# REF# ERR 
I FC  FJ) 25#  20  7 
WRITE!  1#  26) 

FORMAT!  1  3H  ENTER  IDREF) 

READ!  1#*)  IDREF 

GO  TO  ! 30# AO# AO)# IDREF 


WRITE!  1#  32) 

FORMAT! 19H  ENTER  AR  FFR  ALRO) 
READ!  1  #  *  )  AR#  FFR#  ALRO 

PHR=ALRO*Pt/lRO. 

PRI=2«*PI*! 1«+FFR)/SPC 


GO  TO  9  7 
WRITE! 1#  A2) 

FORMAT!  28H  ENTER  BWF1  BWF2  HOLD  CDOFF) 

READ!  l..*)BWFl#BWF2#  HOLD#  CDOFF 

TAU1=1./BWF1 

TAU2= 1 • /BWF2 

BT1=PI*BWF1/!2.*SPC) 

BT2=PI *  BWF2/! 2.*SPC) 

Bl!  1 )  =  <  A**  TAU1-  ?.•  )  *BT1 
Bl  C  2)  =  !  A*  *TAU1“  I  •  )*  BTl 
Bl!3)=!A»ttTAUl+l» )*BT1 
Bl! A)=! A.*TAUl+2.)*BTl 
B2! 1 ) s ! A*  *  TAU2-2# )*  BT2 
B2! 2)=! A.*TAU2-1 • )*BT2 
B2( 3)=! 4.*TAU2+1 • )*BT2 
B2! A)=! A.*TAU2+2.)»BT2 
EX1=EXP!-2.»8T1) 

EX  2=  EXP ! - 2 • *  3T2 ) 

DO  AA  I  =  1 #  A 

EXCOl! I )  =  EXP! “BTl )*COS! Bl !  I  ) ) 

EXC02!I)  =  EXP!“BT2)*C0S!B2!I)) 

0 1 T! I )=0. 

01TT! I )  =  0  • 

02T! I >=0. 

02TT! I >=0. 


DINT=0. 

VI MT=0 • 

JHOLD=HOLD»  SPC+  •  5 


JHL  D=  0 

JDLAY= ! CDOFF* SPC/CODR) + • 
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IFCIDREF-2)97,  50*60 
50  WRITE?  1,52) 

52  FORMAT? I 2H  EM TER  GREF) 

READ? 1,  *) GREF 

SLOPE=GREF»PI*PI*BWFi*BWF2/?0.36«SPC*SPC) 

GO  TO  9  7 

60  WRITE? 1, 62) 

62  FORMAT?  11H  EMTER  ARL) 

READ? 1 ,  * ) ARL 

GLIM=ARL*PI*PI*BWF2/?2.4*SPC) 

9  7  F.J=0. 

CALL  TTYLF?  5) 

CALL  PHS?  1,0) 

GO  TO  ?95>  104,  1  10),  I DREF 
98  CALL  PHS?  2,  2) 

WRITE?  1,  100)1  DREF#  AR,  FFR,  ALRO 

100  FORMAT?  6HI DREF=, I  2*  8X, ‘ 5HFIXED  REFERENCE, /,  3HAR=,F8.4,  5X, 
C4HFFR=,  El 0.4,  2X,  5HALR0=,  F7. 3) 

GO  TO  118 

104  CALL  PHS? 2, 5) 

WRITE? 1, 106) I DREF, GREF 

106  FORMAT?  6HI DREF=, I  2,  8X, 1 6HLI NEAR  REFEREMCE, /, 5HGREF=, F8. 4) 
GO  TO  114 

110  CALL  PHS?  2,  5) 

WRITE?  1,  1  12)  I  DREF,  ARL 

112  FORMAT?  6HI DREF-, I  2, 8X,  1 7HLIMI TED  REFEREMCE, /, 4HARL=, F8 • 4) 
114  WRITE?  1,  1  1  6)  BWF1 ,  BWF2,  CDOFF,  .JDLAY,  HOLD,  JHOLD 

1  16  FORMAT?  5HBWF1=, El  0 *  4, 1 7X,  5HBWF2=» El  0 • 4, /,  6HCD0FF=, F8 »  4, 

C2X,  6H JDL4Y= *  I  6,  /,  5HH0LD=,  F8»  4,  3X, 6HJH0L D=, I  6) 

1 18  CALL  PHS?  3, 10) 

GO  TO  315 

207  GO  TO  ?310, 210,210), I  DREF 
210  DIN=CODLO»SUM 

DO  220  1  =  1,4 

220  01  ?  I  )  =  F.XC01?I)*?2»*01T?  I)-DINT)-EXl*01TT?I)  +  DIN 

OUT1=0. 5»?01? 1) -01? 4)) -01? 2) +01? 3) 

DIMT=DIN 
DO  230  1=  1>  4 
01TT? I >=01T? I ) 

230  01T?I)=01?I) 

GO  TO  ?  310, 240, 250),  I DREF 
240  VIN=SL0PE*0UT1 

GO  TO  280 

250  UIN=  SI GM?  GLIM, OUT  1 ) 

280  DO  290  1=1,4 

290  02? I )=EXC02? I )*?2.*02T? I )- VINT)-EX2»02TT? I )  +  VIN 

OUT2=0 • 5*  ?  02? 1) -02? 4)) -02? 2) +02? 3) 

UINT=VINJ 
DO  300  1  =  1,4 
02TT? I ) =02T? I ) 

300  02T? I ) =02? I ) 
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REF=-CO  JL'J*QUT2 
ERR=REF- SUM 
I F(  JKLD- J40LDD  30  2>  315 
30 R  ERR=0 • 

GO  TO  315- 

310  REF=ARC* SIMCPHR) 
PHR=PHR+PRI 

I  F< PHR-6*  2R31853)  31  3>  31 
3  1  2  PHR=PHR-6«  2P.  31853 

313  ERR=REF-SfP*4 
315  RETURN! 

END 

ENDS 


APPENDIX  III 


In  this  appendix,  the  steady-state  amplitude  of  the  array  output 
signal  which  minimizes  the  mean-square  error  when  the  reference  signal  is 
inverted  for  a  r'raction  of  the  time  is  calculated.  The  steady-state  array 
output  signal  is  assumed  to  be  of  the  form 


(49)  S{t)  =  q  M(t)  sin  Ut  +  $) 

V 

where  constant  q  is  to  be  calculated  and  M(t)  represents  a  +1,  -1 
amplitude  modulation  which  is  equivalent  to  the  actual  bi -phase  (0,180°) 
modulation.  T he  reference  signal  is  assumed  to  have  unity  amplitude  and 
(+1,-1)  amplitude  coding  M(t): 


(50)  R(t)  =  M(t)  sin  (U(.t  +  $). 

The  error  and  squared-error  are  therefore  given  by 

(51)  E(t)  =  [M(t)  -  q  M(t)]  sin  (W(.t  +  +) 
E2(t)  *  [M(t)  -  q  M(t)]2  sin2(u  t  +  <fr). 

w 


The  sine-squared  factor  has  a  constant  time-average  value  of  one-half 
independent  of  q.  To  proceed  further,  the  modulations  R(t)  and  M(t)  will 
be  assumed  to  differ  in  sign  a  fraction  6  of  the  time.  Figure  33  illustrates 
one  possible  distribution  of  the  reference  signal  inversion  periods  con¬ 
tributing  to  the  fraction  6;  the  distribution  does  not  affect  the  average 
square  error  which  is  given  by 


(52)  E2  =  -  {[-(1  +  q) ]2<5  f  (1-q)2  (1-5)} 

=  ~  (4  q5  +  (1-q)2]  . 


The  value  of  q  which  minimizes  the  mean-square  error  is  found  from  the 
equation 

(53)  0  =  1,(46  -  2(1  -,mse>]  ■■ 

<w  * 1  - 25- 


75 


The  application  of  this  result  is  illustrated  by  the  following  three 
examples. 

Example  1.  Square -wave  Data  Modulation 

Assume  that  M{t)  is  an  alternating  +1,  -1  (square-wave;  data  sequence 
which  occurs  at  a  rate  fj  and  that  the  data  filter  delay  equals  one-fourth 
the  data  bit  period.  The  delayed  sequence  is  M(t).  The  signs  of  M(t)  and 
M(t )  differ  during  the  first  quarter  of  each  data  bit;  therefore,  6  =  0.25 
and  qmse  =0.5.  The  suppression  is  6  dB. 

Example  2.  Random  Data  Modulation 

Assume  that  M(t)  and  M(t)  are  random  binary  data  sequences  and  that 
the  data  filter  delay  is  the  same  as  in  Example  1:  one-fourth  a  data  bit 
period.  Since  the  rate  of  occurrence  of  the  sign  differences  is  one-half 
the  data  rate,  it  follows  that  6  =  (0.5) (0.25)  =  0.125.  Therefore, 

%>se  =  which  corresponds  to  a  2.5  dB  suppression. 

Example  3.  Random  Code  Timing  Error 

Assume  that  M(t)  and  M(t)  are  random  binary  code  sequences  which 
are  identical  except  for  a  misalignment  of  their  time  bases  by  one-h^lf 
of  a  code  bit  period:  e  =  0.5  f^.  Now,  5  -  (0.5  fed) (0.5  fC(j)  =  1/4  since 
the  rate  at  which  the  reference  signal  is  inverted  equals  one-half  the 
code  rate.  Therefore,  q  =  1/2;  a  6  dB  reduction. 
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q  M  (t ) 


Example  wa 


APPENDIX  IV 

ADAPTIVE  ARRAY  RESPONSE  FOR  A  CW  INPUT  SIGNAL  AND  CM  REFERENCE 


For  the  special  case  of  cw  input  and  reference  signals,  i.e., 
when  signals  y^(t)  and  R(t)  are  expressible  as 


(54)  y.j  (t)  =  As  cos(wct  +  esi)i  i  =  1 »  —  ,2N 


(55)  R(t)  =  Af  cos(wct  +  er) 

the  system  of  Eqs.  (39)  can  be  written  in  matrix  form  as 


(56) 


dW(t) 


dt 


=  -  A  A 
2  Hs  Hr 


cos(er  -  esl) 


Lcos(er-  es2N) 


K  .2 

1  As 


cik 


w(t) 


Vector  [W(t ) ]  is  a  column  vector  composed  of  the  2N  weighting  coefficients 
and  [C . k]  is  a  2N  x  2N  matrix  with  elements  given  by 

(57)  C.k  -  cos(es.  -  esk)  . 

If  the  incident  signal  arrives  broadside  to  the  linear  array  then 
f  0  i -odd 


(53)  es. 


l -even 


and  the  matrix  is  degenerate: 


(59) 


*ik3  a 


10  10  10 
0  10  10  1 


1  0  1 
0  1  0 


0  1  0 
1  0  1 
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Clearly,  from  (56) 


(60) 


"l  =W3 


_  *  *  • 


=  w. 


2N-1 


so  that 


W2  =  W4 


=  W 


2N 


(61) 


w3(t) 

w5(t) 


W]  (t)  -  W] (0)  +  w3(0) 
W1(t)  -  W^O)  +  W5(0) 


W4(t)  =  W2(t) 
W6(t)  =  W2(t) 


w2(0)  +  w4(0) 
w2(0)  +  w6(0) 


etc. 


etc, 


Substitution  of  (61)  into  (56)  gives  the  two  equations  of  interest: 

•  v  v  9  2N-1 

(62)  W,(t)  =  |  AsAr  cos  8r  -  |  A‘  {N  W,(t)-N  W,(0)  +  \  Wk(0)> 

k-idd 


w2(t)  =  \  AsAr  cos(er  +  f-)  -  |A^{N  W2(t)-N  W2(0)  +  “j  W£(0)} 

£**  2 


K  »2, 


2N 


These  equations  have  solutions  given  by 


NKA 


(63)  W^t)  = 


”Ar 

NA_ 


2N-1 


C0S  9r  +  N  J,  V0) 

k-odd 


«,-even 


A  ,  2N-1 

+  Wq  cos  er  ”  N  Wk(0)  +  Wl(0) 


k-odd 


and 
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(64)  W2(t)  = 


e 


t 


cos(er+  | 


2N 


I 

£=  2 


w*(o) 


£-even 


NKA2 

s 

2 


Ar  w  i  2N 

+  WT  cos(0r+  2  ^  "  N  l  9  V0)  +  W2^‘ 

s  £=2 

a-even 

Equations  (61),  (63),  and  (64)  represent  the  complete  solution  when  the 
signal  arrives  broadside  to  the  array. 


Difficulty  is  encountered  in  attempting  to  obtain  solutions  applicable 
for  arbitrary  angie  of  incidence  due  to  dependencies  among  the  weighting 
coefficient  derivatives.  However,  solutions  have  been  obtained  for  off- 
broadside  incidence  at  angles  corresponding  to  45°  phase  shift/element 
and  90°  phase  shift/element.  The  matrix  Cj|<  is  different  from  (59)  in 
each  of  the  two  cases,  but  it  is  again  degenerate  of  rank  two.  Relationships 
between  the  weighting  coefficient  derivatives  can  again  be  found  and  system 
(56)  reduced  to  two  equations.  Their  exponential  solutions  contain  the 
same  time  constant  as  in  (63)  and  (64)  above.  It  would  appear  that  the  time 
constant  is  independent  of  the  angle-of-incidence  from  these  results. 
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